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Abstract 

In this paper, we investigate the uniform regularity of solutions to the 3-dimensional isentropic 
compressible Navier-Stokes system with free surfaces and study the corresponding asymptotic limits of 
such solutions to that of the compressible Euler system for vanishing viscosity and surface tension. It 
is shown that there exists an unique strong solution to the free boundary problem for the compressible 
Navier-Stokes system in a finite time interval which is independent of the viscosity and the surface 
tension. The solution is uniform bounded both in W^’°° and a conormal Sobolev space. It is also 
shown that the boundary layer for the density is weaker than the one for the velocity field. Based on 
such uniform estimates, the asymptotic limits to the free boundary problem for the ideal compressible 
Euler system with or without surface tension as both the viscosity and the surface tension tend to 
zero, are established by a strong convergence argument. 
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1 Introduction 


Consider the motion of a viscous isentropic compressible fluid with surface tension effect along a free 
boundary. It is governed by the following isentropic compressible Navier-Stokes system written in Eulerian 
coordinate as 


pI + div(p®M®) = 0, 

{p‘^u^)t + div{p'^u^ 0 u^) + = peAu^ + {p + A)eVdivu^, 


( 1 . 1 ) 


where fit is a simply connected domain in occupied by the fluid at time t > 0, p^,u^, which are 
unknowns, represent the density and the velocity field respectively. The pressure is given by the y-law 


7 > 1 


( 1 . 2 ) 


The constant viscosity coefficients pe, Xs satisfy the physical restrictions 


p > 0, ‘2p 3A > 0, 


(1.3) 


where the parameter e > 0 is the inverse of the Reynolds number. We assume the boundary of fit is 
given by 

St = {a; e R3|F'^(a;,t) = 0}, (1.4) 

where F^{x, t) is an unknown function which will be uniquely determined by the velocity field. To study 
the well-posedness theory of this free boundary problem, the following two boundary conditions are 
imposed on St, t > 0. On the one hand, the kinetic boundary condition, which states the fluid particles 
do not cross the free boundary, reads 

dtF^ + ■ V^F^ = 0, for x € Ef (1-5) 


On the other hand, we also need the dynamic boundary condition to balance the stress tensor on the both 
side of the free boundary. When the surface tension is taken into consideration, this boundary condition 
can be written as 

= (2/ieS'u^ + A£divu^)n^ + x € E* (1.6) 

which describes the stress tensor of the fluid is proportional to the mean curvature of the free boundary 
Ej. Here Su^ = + V‘u®) and n® = denotes the outward normal vector of fit, Pe is a given 

constant external pressure, a denotes the surface tension coefficient, F[ is the double mean curvature of 
Ef which can be expressed in the form 

= Asj(t)x, X = (x^,x^,x^), (1.7) 

where {t) is the Laplace-Beltrami operator on Et. We also impose the initial data for the compressible 
Navier-Stokes equations (HH) as 

(p^-«^F^)(0,x) = (pg,itg,T;f)(x), xSHo (1.8) 

such that 

0 < -^ < Po < C'o < oo, (1.9) 

where fig is a given initial domain determined by F^ and Cg is a given constant independent of e. 

The study of fluid motions with free surfaces is an important topic in fluid dynamics. For incom¬ 
pressible viscous fluids, we refer to the works by Beale [J], Tani [iS], Solonnikov [33] and Guo-Tice [3D] 
and references therein for the local well-posedness with or without surface tension. As to compressible 
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viscous fluids, the local-wellposedness theory is established by Secchi-Valli [40] without surface tension, 
and by Solonnikov-Tani [44], Zajaczkowski [^ [^ and Tanaka-Tani [47] with surface tension and ref¬ 
erences therein. For the inviscid fluid, it is much more difficult to get the regularity of free boundary. 
Wu [^ [55] made a big breakthrough for the local well-posedness of irrotational incompressible Euler 
system in two and three dimensional, we also refer to [5] dH [MIIMIISS] and the references therein for 
related works. Later, Wu M proved an almost global existence result in the 2-D case and the global 
existence results in 3-D case is proved by Wu m and Germain-Masmoudi-Shatah[T7| [18]. We also refer 
to [m [m [55] dU [55] and the references therein. For the compressible Euler system, Lindblad [55] proved 
a local well-posedness result by using Lagrangian coordinates and Nash-Moser construction. Using the 
theory of symmetric hyperbolic systems, Trakhinin m provided a different proof for the local existence 
of solutions. Both the estimates in m and m had derivatives loss. Recently, Coutand-Hole-Shkoller 
[13] proved the well-posedness for the motion of a compressible liquid with or without surface tension, 
and with no derivative loss. The zero surface tension limit is also established in [13] . 

Another classical and interesting problem in the mathematical theory of fluid mechanics is to study the 
asymptotic limit of the solutions to the Navier-Stokes equation at high Reynold number which corresponds 
to small viscosity. There has lots of literature on this problem when the domain has no boundaries, see for 
instances [injinKiiEO]. However, in the presence of physical boundaries, the problems become much more 
complicated. When no-slip boundary condition is imposed on the incompressible fluid in a fixed domain, 
the vanishing viscosity limit of the incompressible Navier-Stokes is one of the major open problems due 
to the possible appearance of boundary layers, as illustrated by Prandtl’s theory. In [371 [3H], the authors 
proved the (local in time) convergence of the incompressible Navier-Stokes flows to the Euler flows outside 
the boundary layer and to the prandtl flows in the boundary layer at the inviscid limit for the analytic 
initial data. Recently, Y. Maekawa [27] proved this limit when the initial vorticity is located away from 
the boundary in 2-D half plane. While, for the incompressible Navier-Stokes system with Navier-slip 
boundary condition in a fixed domain, considerable progress has been made on this problem. Indeed, the 
uniform bound and a uniform existence time interval as e tends to zero are obtained by Xiao-Xin in 
m for flat boundaries, which are generalized to in mm- However, such results can not be expected 
for general curved boundaries since boundary layer may appear due to non-trivial curvature as pointed 
out in |21j . In such a case, Iftimie and Sueur have proved the convergence of the viscous solutions to the 
inviscid Euler solutions in L'^(0, T,; T^)-space by a careful construction of boundary layer expansions 
and energy estimates. However, to identify precisely the asymptotic structure and get the convergence in 
stronger norms such as L°°{0,T-, H^){s > 0), further a priori estimates and analysis are needed. Recently, 
Masmoudi-Rousset m established conormal uniform estimates for 3-D general smooth domains with the 
Naiver-slip boundary condition, which, in particular, implies the uniform boundedness of the normal first 
order derivatives of the velocity field. This allows the authors to obtain the convergence of the viscous 
solutions to the inviscid ones by a compact argument. Based on the uniform estimates in [31] . better 
convergence with rates have been studied in m and [55]. In particular, Xiao-Xin [55] has proved the 
convergence in L°°{0,T-, H^) with an rate of convergence. 

For the compressible Navier-Stokes equations in fixed domain, Xin-Yanagisawa |61j studied the van¬ 
ishing viscosity limit of the linearized compressible Navier-Stokes system with the no-slip boundary 
condition in the 2-D half plane. Recently, Wang-Williams m constructed a boundary layer solution of 
the compressible Navier-Stokes equations with Navier-slip boundary conditions in 2-D half plane. The 
layers constructed in m are of width 0{^/e) as the Prandtl boundary layer, but are of amplitude 0{^/e) 
which is similar to the one [21] for the incompressible case. So, in general, it is impossible to obtain 
the or W^’^{p > 3) estimates for the compressible Navier-Stokes system (11.111 with the generalized 
Navier-slip boundary condition. Later, Paddick [35] obtained an existence and conormal Sobolev regu¬ 
larity of strong solutions to the 3-D compressible isentropic Navier-Stokes system on the half-space with 
a Navier boundary condition. Recently, Wang-Xin-Yong [53] also obtained an uniform regularity for the 
solutions of the compressible Navier-Stokes with general Navier-slip boundary conditions in 3-D domains 
with curvature, especially, the vanishing viscosity limit of viscous solution to the corresponding inviscid 
one was also obtained with rate of convergence in L°°. In [53], it is also shown that the boundary layer 
for density is weaker than the one for velocity fields. 

As to the vanishing viscosity limit problem of fluid motion with free surfaces, an interesting result was 
achieved recently, by Masmoudi-Rousset [35], that the local existence of solutions to the incompressible 
Navier-Stokes system with the gravity field but without surface tension in an uniform in e time interval 
by using a suitable functional framework based on conormal Sobolev spaces which minimizes the needed 
amount of normal regularity but which gives a control of the Lipschitz norm of the solution. These 
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regularities of solution are necessary and reasonable for the free boundary problem. Actually, the Lipschitz 
regularity guarantee that the free surface moves along the partical path into outside of the fluid. And, 
in the vicinity of the free boundary, the expected behavior of the solution to the incompressible Navier- 
Stokes equation is u^(x^ t) ~ u{x, t) + {t, y, -^) so that it is hard to obtain an uniform iJ^-norm(fc > 

2) estimate for in a time interval independent of e. In particular, a new existence result for the 
incompressible Euler equations can be obtained by strong convergence argument for the Navier-Stokes 
equations. Later, Elgindi-Lee has tried to study the similar problem in the presence of the surface tension 
in m- It should be noted that in the case surface tension be a positive constant, the pressure term in the 
Euler system becomes less regular. On the other hand, it is also interesting to investigate the zero surface 
tension limit for free boundary problems. This is nontrivial since that the surface tension coefficient a is 
connected to the mean curvature of free surface which is a second order derivative term of the boundary 
function. It is a subtle issue to obtain the uniform estimate in a. Ambrose-Masmoudi ia[3] studied the 
zero surface tension limit of irrotational water waves. The zero surface tension limit of incompressible 
Navier-Stokes equation has been established by Tan-Wang in for small initial data. For the zero 
surface tension limit of compressible Euler with free surfaces, we can refer Coutand-Hole-Shkoller |13j . 

In this paper, we are interested in the existence of strong solution to the compressible Navier-Stokes 
system dm with free surfaces in a finite interval of time independent of viscosity e and the surface 
tension a and the corresponding asymptotic limits problem of vanishing viscosity and zero surface tension. 
Formally, when the £ tends to zero, the limit of (HH) is compressible Euler system with surface tension 


dtp + div(pM) = 0, 

pdtu + p{u ■ V)u -1- Vp = 0, 


(I.IO) 


with the boundary conditions 

dtF + u-\7F = 0, and p = Pe — (TH, x£T,t, (l-H) 

Moreover, when the surface tension coefficient tr also goes to zero, it can be reduce to the compressible 
Euler system (11.1011 with the boundary conditions 

dtF + u-\7F = 0, and p = pe, a; G S*, (1-12) 

Motivated by Masmoudi-Rosset [32] . we aim to obtain the uniform in e and cr regularity in the anisotropic 
conormal Sobolev spaces and a control of the Lipschitz norm for solutions to the compressible Navier- 
Stokes equations (EU with free surfaces. Based on this uniform regularity, we can study the asymptotic 
limits to the ideal compressible Euler system with free surfaces as both the viscosity and the surface 
tension tend to zero by a strong compactness argument. Of course, the Taylor sign condition for the 
pressure on the boundary is needed to have local well-posedness for the Euler system. Usually, one 
expects that the boundary layer for density is weaker than the one for the velocity fields, we also aim to 
prove this facts in this paper. It should be mentioned that around same time as this research project, 
similar theory has been just established by Wang-Xin[5^ for the incompressible fluids with rigid bottom 
below. Although there are some common difficulties caused by the appearance of surface tension for both 
cases, yet, as we shall explain later, there are major differences in these two cases. Different techniques 
are needed in the compressible case due to the lack of divergence free conditions which are essential in 

[221 [Sg. 


1.1 Reformulation the problem into local coordinates 

Due to the possible appearance of boundary layers, we shall consider this problem in the conormal Sobolev 
space as [32j . In order to define the conormal derivatives for the free boundary problem, we hrst assume 
that the initial domain Hq has a covering such that 

Uq (e Uq ^k=i (1-13) 

with each fc = 0, • • • ,n being a convex domain. Here, Dq Hq is chosen to satisfy 


dist(Eo,Ho) = 2(io > 0. 


( 1 . 14 ) 
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where Sq is the boundary of flo and do is a given positive constant. For each fixed k and a > 0, define 
^k,a as 

^fc,o = {x G I dist(x, fjfc) < ado}. (1-15) 

Clearly, 

ilo <E i^o,i fifcp. (1-16) 

^k,a,k = 1, • • • ,n will be called the boundary covering. Without loss of generality, we assume that 

rio n Ctk,i = {x|x3 < hi{xi,x 2 )} n nk,i, dfto n nk,i = {x|x3 = hi{xi,x 2 )} n 


since the other cases can be handled similarly. Define 

= {(xi, X2) I X G n Do, for some x = (xi, X2, X3) G M^}, fc = l,---,n. 


Then, is the projection of ilk,a H Dq onto the hyperplane x R^j^. 

Similarly, for Dj, we can define h\{t, xi,X 2 ) and E( ’“. 

Remark 1.1 Under the assumption of the velocity is bounded, it is easy to know that Dpj, D,;. 1 , k = 
1, • • • ,n is still an effective covering of Dt at least in a short time t G [0, Tg], i.e. 


Dt (s Do.i Dfe i, Vt G [0,ro], with Tq < 1, 
where Tq depends on the upper bound of velocity and do ■ 


As in m, to deal with the free boundary problem, one can reduce the problem in each boundary 
covering into a fixed domain and use the local coordinate to define the conormal derivatives. Instead of 
using Lagrangian variables, we define a family of differmorphism $|(t, •) by the following process. Since 
hf, is defined locally, we first extend hf, to R^ by multiplying a cut-off function 


hUt,y) = fik{y)hl{t,y), 
where ifkiy) is a smooth function satisfying 


Aiy) = 


1; 2/ £ (0), 

0, 2/ G 


Then, define as 

= K{z{^))hl{t,f), 

where " denotes the Fourier transform with respect to y and k is defined by 


(1.17) 


(1.18) 


(1.19) 

( 1 . 20 ) 


which implies that 

yk(t,y^z) = ^e-^^ f ■^e^^^h%{t,y)dy. ( 1 . 21 ) 

z J-g2 Z 

As in [32], (/j| can be chosen as 

^Ut:y,z) = Az + ylit,y,z), (1.22) 

where A > 0 is a constant chosen such that 


dz<fl{0,y,z)>l, {y,z) eS. (1.23) 

Therefore, the differmorphism <1>| can be defined by 

$|(t, •) : 5 = R^ X (- 00 ,0) Dk{t) D Dt n 

{y,z)^x:={y,ipl{t,y,z)), (1.24) 


where Dk{t) is defined by 


Dk{t) := {x = {y,X 3 ) | j/ G R^ X 3 < h%{t,y)}. 


(1.25) 
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By the above parameterizations, the free surface becomes z = 0 locally. Although we have defined an 
extension of h% in x (—00,0), but it will be more convenient to work in a vicinity of the boundary, i.e. 


Sk,aii) = {(y> liVy S such that x = (y, y, z)) € C 5, fc = 1, • • • , n. 

with 0 < a < Since the velocity will be uniform bounded, one has that 


(1.26) 


C 5,,1(0) C C 5,,1(0) C C 5,,|(0), Vt e [0,ro], To « 1. (1.27) 

Since Ej is given locally by X 3 = h{t, y) (Henceforth the subscript k and the supscript e in p^,u^, h^, 
will be omitted for notational convenience), it is convenient to use the local coordinates (y, z) defined in 
(11.241) at least locally. Once the choice of (p is made, based on the observation of (11.271) . one can reduce 
locally(neighborhood near the boundary Tit) the problem into the fixed domain 51(0) by setting 

0 {t,y,z) = p{t,<^{t,y,z)), v{t,y,z) =u{t,<^{t,y,z)), {y,z)eSi{ 0 ). (1.28) 

For simplicity, we will denote also {p,u) as {g,v) in the interior domain. 

Now, we rewrite the system dm) to the one for (g,v,h) in the new coordinate {y,z) of the fixed 
domain 5i (0). First, one introduces the operators , i = t, 1, 2, 3 such that 


dtQ= {d,p) o dfv = (5*m) o (t,$(t, •)), 


which directly yields that 


= * = t,i,2, 

,^3 = 


(1.29) 

(1.30) 


Then, the compressible Navier-Stokes system dm is locally reduced to the system for (p, v) as 


df g + div‘^(pv) = 0, 
gdfv + gv ■ V'^v + V'^p{g) 


2pedvv^{S'^v) + XeV'^div'^v, 


{y,z) € 5i(0), t > 0, 


(1.31) 


where div‘f’v = ^2 > v-V’^ = and S'^v = \{y^v + V^v^). And the 

two boundary conditions read as 


dth = v'H = -Vy{t, y, 0) • \7yh + v^it, y, 0), {y, z) e 5i (0) (1.32) 

and 

pN = (2^£5‘^^ + A£div‘^z;)N + p^N -aVy , (y, 0) e 5i (0), (1.33) 

\/l+|Vy/l|2 

where v = {vy,V 3 )* := (rii,U 2 ,^ 3 )*, V = {Vy,dzY '■= {di,d 2 ,dzy and N = {—dih,—d 2 h,\y. For later 
use, we dehne n = . 

Remark 1.2 By using the covering of fit, one can always assume that each vector field {g,v,h) is 
supported in either 5^. i(0) in {y,z) coordinate, or in Hop, since this assumption can be achieved by 
multiplying a cut-off function ipj. € (^((“(Mi).) such that ipj. = 1, iy,z) € 5;. i(0) and ifk = 0, {y,z) € 
— 5^. 1 (0), and ipo G C((°(M^) such that tpo{x) = 1 for x G Hq.I f’oix) = 0 for a: G — Hop, 
respectively. Therefore, using this localization arguments, one can also assume that the problem p.31|) 
holds in S or and (ll.32|) - (ll.33l) hold on z = 0, ?/ G M^. 

To measure the regularity of functions defined in S, we shall use the Sobolev conormal spaces as in 
[32] . Introduce the vector fields 

Z, = d,,i = l,2, Z3 = -^dz. (1.34) 

i — Z 

The Sobolev conormal spaces is defined as 

= {/ G L\S) I Z“/ G L^S), |a| < m}. 


(1.35) 
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where with norm defined as 

ll/ll™ = E ii^“/ii'- 

|Q;|<m 


Similarly, set 

W^,^{S) = {f &L°°{S)\Z‘^f &L°°{S), |a|<m}, (1.36) 

and 

ii/ikoo= E ii^“/iU“- 
|Q:|<m 

However, for the compressible Navier-Stokes, the time derivative should be involved in the energy 
functional. So, we also set 

Zo = 5o = dt, (1.37) 

and, for smooth space-time function f{t,y,z) with t > 0, (y^z) € S, set the following notations 


Wfmnr. = E ik“/(of- ii/wii«.,=o = E ik“/(t)iil 

|a|<m 


||VZ™/(t)f = ^ ||VZ“/(t)f, 

|Q;|<m 

and 

i/(i)ik= E k“/(t,-,^ = o)k, \2^fit)\i= E k“/(^,•,^ = o)|^ 

|Q!|<m |a|<m 


Remark 1.3 It should emphasized that the norms defined above can be used in each covering although 
we have not written this fact explicitly. 

Remark 1.4 In the present paper, we will focus on the estimates near the boundary, i.e. in iSi(O). The 
interior estimate is easy to obtain in a similar way by a cut-off approach such that no boundary terms 
involves, since the conormal Sobolev space is equivalent to the standard Sobolev space in the interior 
domain. The global estimates follow from collecting all the local boundary ones and the interior one. 

Notations: Throughout this paper, the positive generic constants that are independent of £,a are 
denoted by c, C. And || • || denotes the standard norm, and || • ||_ym (to = 1,2,3,-••) to denote the 
Sobolev iJ™ norm. The notation | • |s will be used for the standard Sobolev iJ® norm of functions defined 
on boundary, and note that this norm involves only tangential derivatives, a < b used here denotes 
a < Cb for some positive constant independent of e,cr. A(-) denotes a polynomial function independent 
of £ and cr, and may change from line to line. 


1.2 Main results 

The aim of this paper is to get a local well-posedness result for strong solutions to (ED in an interval of 
time independent of e and cr for eg (0,1], tr S [0,1]. Note that such a result will also imply the local 
existence of strong solutions for the Euler equations with or without surface tension. As it is well-known 
[H ESI SUES], in the absence of viscosity and surface tension, we need to require the initial pressure po 
for the compressible Euler flow satisfies the Taylor sign condition (see, for example m and m), i-e- 

dj)^ 

0<co<-^, x€Et. (1.38) 

In the local coordinate (y, z), the initial data (11.811 is reduced to 

(^'^ z, 0) = (yg, /i§)(y, z), (1.39) 

0 < 7 ^ < pg < Co < oo. 


such that 


(1.40) 


and the Taylor sign condition becomes 


0 < Co < -O'^Pq, z = 0. 


(1.41) 


Since boundary layers may appear in the vicinity of physical boundaries, in order to obtain uniform 
estimates for solutions to the compressible Navier-Stokes system with free surfaces, one needs to hnd a 
suitable functional space. Here, we define the functional space X^{T) for (p, v, h) = {p, v, h){t) as follows: 


X^(T) = ^{p,v,h) G L°°{[0,T],L'^); esssupo<t<Tll(p,'c, ^)(i)||x^ <+oo|, (1.42) 

where the norm IK-, •, is given by 

\\iP,v,h){t)\\l-^ = \\ip,v){t)\\l^,n + \\V{p,v){t)\\'^m-2 + ||Ap(t)||^i + |(h,v^Vy/i)(t)|^„ 

+ II V(p, vmWn^,^ + ell V(p,u)(t)||?,„_. + e|| Ap(t)||?,. + e|| V\(t)||i,. + e|Z-h||. ( 1 . 43 ) 


Note that the a priori estimates in Theorem 14.11 below is obtained in the case that the approximate 
solution is sufficient smooth up to the boundary. In order to obtain a selfcontained result, one needs to 
assume that the approximate initial data satisfies the boundary compatibility conditions. Therefore, we 
set 


-^NS,ap 


|(p, m) S C'^"‘(Ho), T'(x(si, S 2)) G The Taylor sign condition (11.381) holds; 


dfU and d^F, k = 1, • • • , m are defined through the Navier-Stokes equations 
(HU and (HU, respectively; d^{p, u, F), fc = 0,1, • • • , m — 1 satisfy the 

boundary compatibility conditions (11.61) >, (1-44) 


where the domain Hq is defined through F as previous. Note that the definition of is given in 

the Euler coordinate, but it will be changed automatically into local coordinate when one evaluates it by 
some norm, and we should keep this in mind all the time. Define 


A 


£,m _ closure of ^^5 in the norm IK-, •, • 


NS 




We assume that the initial data (pg,ug, hg) G X'^^ and satisfies 

sup {\\iPo,Vo,hl)fxe + ||Vz;g||^™_i| < Co, 

£6(0,1], <tG[0,1] 


(1.45) 


(1.46) 


where pg = p(Po)i Co > 0, Cg > 0 are positive constants independent of £ G (0,1] and a G [0,1]. Thus, 
the initial data (pg, Ug) is assumed to have higher space regularity and compatibilities. It should be noted 
that the initial data may also depends on tr, but, for simplicity, this will not be written explicitly. 

Our main result is as follows: 


Theorem 1.5 (Uniform Regularity) Assume that m > 6, e G (0,1] and a G [0,1]. Consider the 
initial data (pg,Mg,Eg) G given in (11.81) satisfying (11.401) . (11.411) and (ll.46|) . Then there exists a 

time Tq > 0 and Ci > 0 independent of e and a, such that there exists an unique solution {p‘^,u^, F^^) 
to the free surface problem (11.11) . (11.51) and (11.61) which is defined on [0,To] and satisfies the following 
estimates: 


sup ||(p^u^)(^)||^. + ||V(p^u^)(^)||^._. + ||Ap-(t)||^. + 

0 <t<To I 

+ ||V(p^«^)(^)||?,co. +£||V(p^U^)(t)||?,_. +£||Ap^(t)||?,. +£||VV|li,. +£|Z-h|| 


rTo 


rTo 


+ / II Vp^(t)||^.-. + II Ap^(t)||^.dt + / \\Xv^mnr.-.dt 


rTo 


rTo 


+ e ||Vu^(t)||^.+ ||VV(t)||^_.dt + £M ||VV(t)||^_iA<Ci, 


(1.47) 


and 


T^<g%t)<2Co VtG[0,To], 

ZLyQ 


(1.48) 


where Ci only depends on Co, Co and do. 
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Note that in the above result, time derivatives shall be involved in the energy functional. This is 
not necessary for incompressible flow without surface tension studied by Masmoudi-Rousset[32], since 
the pressure p there solves an elliptic equation with a Dirichlet boundary condition so that elliptic 
regularity of p in the spatial conormal spaces can be used to obtain the uniform estimates. But, as to 
the incompressible flow with surface tension studied by Wang-Xin [S2], time derivatives are also needed 
since the Neumann boundary condition for p, which was derived from the momentum equation, shall be 
used. It should be noted that the main difficulties in |52j for the incompressible flow with surface tension 
are that there are no any estimates of the highest order time derivatives d^^p. However, the difficulties 
for the compressible fluid are different, since the pressure p here satisfies a transport equation so that the 
regularity of can be derived by an energy method. As will be seen later, the main difficulties here 
are to derive delicate estimates of div'^n and uniform in a estimates (when cr ^ 0) of the boundary term 
as (14.2.151) involving surface tension. We also give some remarks on our uniform regularity in Theorem 
1.5 as follows. 

Remark 1.6 Note that {p^, F^) in Theorem 11.51 is the solution to (11.11) (11.51) and (11.61) given in 
Eulerian coordinate of Tit- While, {g'^,v‘^, h'^) defined locally in the fixed domain S is used in all the a 
priori estimates. Since the solution has enough regularity, they can he changed to each other, eguivalently. 
These two notations will be used without confusion throughout this paper. 

Remark 1.7 The Taylor sign condition (11.381) or (11.411) is necessary when one studies both vanishing 
viscosity and zero surface tension limits to the ideal compressible Euler system. However, when surface 
tension is considered(for fixed a > 0),this is not necessary since the surface tension prevent the Rayleigh- 
Taylor instability 

Remark 1.8 One can construct a class of initial data to satisfy the compatibility conditions (II .44|) . e.g. 
Let Ho = R_Ro be a ball with radius Rq. When a > 0 is fixed, let (pg, wg) be sufficiently smooth functions, 
and pg be a positive constant poc and ng vanishes in a vicinity of the boundary. By choosing p(poc) = 
Pe + then it is obvious that (pg,ng) satisfies the boundary compatibility conditions. When a = 0, 
it becomes a little complicated due to the Taylor sign condition. In this case, one can assume that the 
initial data is radial symmetric (near the boundary of ball), which are polynomial functions of radius 
Rq — r, so that the boundary compatibility conditions can be reduce to algebraic eguations for coefficients 
of polynomials. We omit the details here for simplicity. 

Remark 1.9 It follows from the uniform estimate of \\Ap‘^(t)\\'^i that the boundary layer for the density 
p^ is weaker than the one for the velocity u® as expected. 

Based on the uniform estimates given in Theorem 1 1.51 we can justify the vanishing viscosity limit, zero 
surface tension limit and the local existence of solution to the free surface problem for the compressible 
Euler system. 

Theorem 1.10 (Inviscid Limit) For any fixed a > 0, under the assumptions in Theorem \l.fA and in 
addition that there exists {go,uo,ho) such that 

lim |K£'o - 00,'^0 - vo)\\l'^ + |hg - /io|l2 = 0. (1.49) 

£—>■0 

Then, there exists {p,u, F)[x,t) on the time interval [0,ro], and satisfies 

sup I ||(p,n)(t)||«m + ||V(p,z;)(t )||^„-2 + \{h, yh){t)\j^r„ + ||Ap(t)||^i + ||V(p, n)(t)||^i,oo 1 

0<t<To I J 

+ f ||Vp(t)|||^m-i(it + f IIAp(t)||^ 2 dt + f ||Vn(t)||^m-idt < Cl < oo, (1.50) 

Jo Jo Jo 

0 < < p(t) < 2Co < oo, Vte[0,ro], (1.51) 

and 

lim sup ||(p^ — p, — i;)(t)||Loo + |h^ — h|iyi,oo = 0, (1-52) 
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where (g,v,h) is the localized representation of {p,u, F). Furthermore, {p,u,F) is the unique solution to 
the free surface Euler equations 


dtp + div{pu) = 0, 

pdtu + p{u ■ V)m + Vp = 0, 


(1.53) 


with the boundary conditions 


dtF + u-VF = 0, and p = Pe — aF[, xGHt, (1-54) 

where fit is the domain occupied by the fluid on time t > 0 and the boundary of Qt is given by 

E,t = {x G F{x,t) = 0}. 

Similarly, by the strong compactness argument in Theorem ll.lOl as £ and a tend to zero independently, 
one can obtain that 

Theorem 1.11 (Inviscid and Zero Surface Tension Limit) Under the assumvtions of Theorem \1.5[ 
if we assume in addition that there exists {go,uo,ho) such that 

lim WiPo^" - go,VQ’‘" -vo)\\l^ + \hl’‘'- ho\L2 =0. (1.55) 

6 — 

Then, there exists {p,u, F)(t) on the time interval [0,ro], and satisfies 

sup ( I|(p,u)(t)||^. + ||V(p,u)(t)||?,_. + \h{t)\l„. + WApmlt + ||V(p,u)(t)||?,.„. j 
0<t<To I J 

+ f ||Vp(t)||^m_i(it + f \\Ap{t)\\'^2dt + [ \\S/v{t)\\'^,n-idt < Cl < oo, (1.56) 

Jo Jo Jo 

0 < :^ < p(t) < 2Co < oo, Vte[0,ro], (1.57) 

ZOo 

and 

lim sup - £», - ■!;)(t)||Loo + -/i|^yi,oo = 0, (1.58) 

®-*-Ote[o.To] 

where {g,v,h) is the localized version of {p,u, F). Furthermore, {p,u,F) is the unique solution to the free 
surface Euler equations 


dtp + div{pu) = 0, 

pdtU + p{u ■ V)u + Vp = 0, 


(1.59) 


with the boundary conditions 


dtF + u-\7F = 0, and p = Pe, x G St, (1.60) 

where fit is the domain occupied by the fluid on time t > 0 and the boundary of Qt is given by 

St = {x G F{x,t) = 0}. 

Remark 1.12 Based on the uniform regularity (ll.47|) - (ll.48l) . for any fixed £ > 0, one can also obtain 
the zero surface tension limit of free surface compressible Navier-Stokes system with surface tension by 
the similar strong compactness argument above. 

1.3 Sketch of the proof 

Since the classical local existence results to smooth solutions of (HU are available [40l |63], the main 
difficulty in the proof of Theorem 11.51 is to get the a priori estimates of the solution in a small time 
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interval independent of e and cr. For notational convenient, we drop the superscript e in the a priori 
estimates. Define 

0m{p,v,h){T) = sup (1 + ||(p,'y,h)(t)|||.e ) + [ ||Vp(t)||^™_i + ||Ap(t )||^2 + ||Vn(t)||^™_idt 
0 <t<T ™ Jo 

+ £ /* \\Vv{t)\\l^rn -\-\\V‘^v{t)\\l^rn-2dt-\-e‘^ f 11 V (i) 11 ^ 1 < +00. 

Jo Jo 

The major goal of this paper will be to derive a uniform bound of the above quantity in a finite time 
interval. We outline the main steps and ideas for the a priori estimates as follows: 


Step 1: Conormal energy estimates for {p,v,h). The first step will be to estimate Z°‘{p,v,h) for 
|a| < m. The basic energy estimate(a = 0) of {p,v, h) is easy to obtain from the total energy identity of 
compressible Navier-Stokes system. In order to get the estimates for higher order conormal derivatives, 
one can perform the Lf energy estimates on the Alinhac good unknowns <5“ = Z°‘p — dfpZ°‘r], = 
Z‘^v — dfvZ°‘ri for |q;| 7 ^ 0 to overcome the loss of ^ derivative for h on the free surfaces as in [52], when 
O’ = 0. However, one subtle difficulty arises on the uniform estimate of boundary terms for and 
This is because the following boundary estimate 

^ A(||u||i,oo + |^|2.oo)(|?^|s+l + |^|s+l), (1-62) 


which plays an key role in the analysis of [32], is invalid for the compressible case due to div'^’u is not 
free. Here and henceforth, ()^ denotes the function taking value on the boundary z = 0. So, in general, 
the term || V^u||^„_i(iT will appear on the right hand side of the estimate for and when one 
applies directly the trace estimates to terms involving f* |(9zf)^|^m-idT(cf. Lemma 1531) . Fortunately, 
with a there, this term can be bounded by the dissipation terms in the estimates for normal derivatives 
and div'^’u in the next steps. When the effect of surface tension((T 7 ^ 0) is taken into consideration, the 
boundary term involving the mean curvature ctV„ • , of the free surface seems difficult to be 


bounded uniformly , since 




Vyh 


V^ + \^vh\' 


= (a™ ■ dtS/yh)dy ^ or a\v' 


..b\2 

\H" 


(1.63) 


which both terms are difficult to be bounded uniformly (cf. Lemma 13. Sp . To deal with this boundary 
term, we will take full advantage of the third component of dynamic boundary condition (13.2611 and 
momentum equations (EU) to turn it into the volume integral (cf. (l4.2.16p - (l4.2.2ip L The key observation 
is that the worst boundary terms appearing in (14.2.1711 and (14.2.2111 by this approach will be canceled 
with each other. It should be reminded that, for any fixed ct > 0, one can derive the uniform in e and a 
estimates without using the Alinhac good unknowns since the boundary regularity can be improved due 
to the surface tension effect. But, when one studies the zero surface tension limit problem, this improved 
regularity of boundary disappear so that the Alinhac good unknowns should be used here. Precisely, one 
can obtain that 


||(H-,Q-)(<)|| + |(/r,v^V,h)(t)||,„.+e / 

^0 

<A{Co) + S f ||Vp||^™_i+£2||v2u||^„_idT + A(0™) [ \\Vv\\l^^-^dT + CstA{0rn)■ 


>0 
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provided the Taylor sign condition (11.4111 holds true. 


Step 2: Estimates for (V^p, div^v). Next, one needs to bound Vp to close the argument. We shall 
use an energy method to estimate Vp, since the pressure p = , shown in (II.31^ ^ for compressible flows, 

satisfies a transport equation due to (II.3111 ^. Thus, one also needs to bound div'^n. This is different from 
the estimate of pressure in [321152] for incompressible flows, which satisfies an elliptic equation due to the 
divergence free condition. Furthermore, since 

dzV = div'^u + n(i92u) + VyU, 
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the estimate of is also helpful to the estimates of normal derivative dzV. It follows from the standard 
energy method that 


|Vu( 




+ e||(div‘^u,V‘^p)(t)||^^„i + f e\\V^div‘^vf+e^\\V^Z^-^dw‘^vfdr 

Jo 


<A(C'o)+feM ||V^u||^„_idT + A(0„) / ||Vu||^„_idT + C 5 tA( 0 „), 


||Vp(t)||^™_.+ / ||Vp||^^_idr< A(Co)+eM ||+ tA(0„), 

Jo Jo 

lldiv'^ull^m-i < A{Om)- 

Here the estimate of ||(Vu, Vp)(t )||^^_2 is a consequence of the fundamental theorem of calculus. It 
should be reminded that the estimate of 

e||div‘^u(t)||^^_i + f e\\V^Z"^-^div^vf + V^Z'^-^iv'^uf dr (1.64) 

Jo 

is the key to bound e||Vu(t)|||^m_i + Jq e\\V‘^v\\‘^m -2 + e‘^\\V^v\\'^rn-idT. This is not necessary for incom¬ 
pressible case, since /q e|| V^u||^m _2 + e^|| V^u||^m_idT does not appear due to (ll.62|) . 


Step 3: Normal derivative estimates Part I. Since VyU and div"^?; have been bounded in the above 
steps, in order to bound d^v, it remains to estimate n(i92;u), where H = Id — n® n denotes the tangential 
vector field. One can follow the arguments in m to derive the estimate for 

Sn = n(5'^uN), 

which is equivalent to Ti{dzv) but vanishes on the boundary. The main conclusion of this step is 

\\Sn{t)f ^^.2 + e\\Sn{t)f^^., + [ e\\V^Z^-^Snf+e^\\V^Z^-^SnfdT 

Jo 

<A(Co) + A(0™) [ \\Vvf^m-^dT + CstA{0m). 

Jo 

This, together with (11.641) . yields immediately that e||Vz;(t)||^,„„i -I- /g e||V^?;||^m _2 + e^||V^u||^,„„idr, 
which does not appear in the arguments of [32] due to the divergence free condition. It should be noted 
that e:||(Vu, Vp)(t)||^m_i is necessary to close the estimate of /p £r^|| V^u||^m-idT. 


Step 4- Normal derivative estimates Part II. In order to close the a priori estimate, it remains to 
bound the m — 1 order conormal derivatives of d^v. As indicated in [32] . one can only expect the estimate 
of /g II Vu||^m_idT. Based on the estimate of div'^u in Step 2 and the fact d^v = div'^u -|- w x N -(- VyV, it 
suffices to bound /g ||a; x N||^,„_idr. Due to the less regularity of Vp and div'^u, one chooses to bound 
U! X N as in [32] instead of S'n, since V^p and V'^div'^u can be eliminated in the equation of the vorticity 
u!. We will follow the approach of microlocal symmetrizer in [32] to derive this kind of estimate. However, 
if one applies directly this argument to f* ||u;||^m_icir as m so that 

(^J ||w||«m-ldT^ <A(Co)+\/ey \u^{T)\l^m-ldT + ti A{Orn)- 


Then it seems difficult to bound y/e f* \uj^{T)\'^rn-idT, since it involves \(dzv)^\'f^,n-i but in general (ll.62p 
is not valid for the compressible flows. It should be noted that even though it follows from trace theorem 
that 


Vs [ |w(t)|^, 

but, eJ^Vwll^™ 
and (I4.5.2|) l that 


pt 1 1 

-idT<^/e \\Vu}\\u^-i\\u}\\u^-idT+ VA{em) < e || Vw||^m-idT-h A( 0 ^), 

Jo Jo 

idr may not be expected to be bounded. However, direct calculations yield(c.f. (|4.4.5p 


(w X N)'' = -2n{{div ■ N, d2V ■ N, 0)‘}. 
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So, one can expect to bound ||a; x N||^m_i(ir. To avoid involving too much regularity of h{due to the 
commutator for x N ), it is equivalent to bound /g \\Z'^~^uj x N||^c?t base on the fact 

[ ||w X N||^™_idT < f \\Z"^-^UJ XNfdT + tA{Om). 

Jo JO 

Note that the boundary condition (11.3311 yields 

X N)*” ^ Z^-^Vyh + {Z'^-^Vyvf + {Z^-^Lof + l.o.t. 


Then, it follows from the argument of the microlocal symmetrizer in |32j that 

[ \\uj X ■N\\^^.^dT < A{Co) + V^AiOm) [ \{Z"^-^uj X ■N)\l,dT + tAiOm) 


< A(Co) + e / llVy'^r + ||V^n||^„_.dr + tA(0^). 


Therefore, the conclusion of this step is 

rt 


< A{Co)+sJ^ ||Vl/™f + ||V\||^_.dr + tA(0„). 


Step 5: L°°-estimates. Finally, we will bound L°“-norm of ||(Vp, Vz;)||^i,oo and i/ell V^w||loc, to close 
the priori estimate. In order to estimate || Vp||^i,cx,, using the anisotropic Sobolev inequality, one only 
needs to estimate ||A‘^p||^i. Actually, one can obtain, for m > 6 , that 


(||AXt)||^, + e\\A^p{t)f^,) + \\Z<^A^pfdT < A(Co) + tA(0™). 

This estimate also implies that the boundary layer for density is weaker than the one for velocity. 

In order to estimate ||Vn||^i,oo and i/e|| , it suffices to to estimate ||5'n||^i.oo and e|| VS'n|||cx,, 
since 

\\Vv\\\^.^ + ellV^lli,. < A(Co) + \\A^P\\\^ + ||^n||?,uo. + e||V^„||i.. 

However, ||«S'n||^i,oo and e|| VS'nllicx, seems difficult to be bounded if one performs L°° estimates directly 
on the convection diffusion equation solved by due to the appearance of the term ||Vdiv‘^ti||^ 2 ,oo, 
which is difficult to controll. Therefore, we try to estimate 

Cn = n(a; X N) - 2T\{{A/y,Qfdtii - (V‘"N)‘u}, 

which also solves a convection diffusion equation but involving || Vdiv''’z;||^i,oo only. By using the similar 
argument in m, one can obtain, for m > 6, that 

+ ellV\||i^ < A(Co) + t^AiOm). 

The estimate of 0m on some uniform time follows by collecting the estimates of the five steps. Note 
that at the end, we need to check that the Taylor sign condition and <I> being a differmorphism remain 
true. 


The rest of the paper is organized as follows: In the next section, we collect some elementary facts and 
inequalities to be used later. In section 3, we study the equations satisfied by {Z°‘p, Z°‘v, Z°‘h) and prove 
the estimates for the lower order commutators and boundary values. In section 4, the a priori estimates 
Theorem I4.1I will be established, which is the main part of this paper. Using the a priori estimates, 
Theorem 1 1. 5I and Theorems II. 101 fl. Ill are proved in section 5. The Appendix collects a generalization of 
Lemma 14 and Lemma 15 for L°° estimates in [3T] so that it can be applied to compressible Navier-Stokes 
equations. 
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2 Preliminaries 


In this section, we will list some basic properties of conormal function spaces and some elementary facts 
which will be used frequently in the process of a priori estimates. For later use, define 

W™([0,T] X 5) = {f{t,y,z) € L^{[0,T] x 5)| Z“/ e L^{[0,T] x 5), |a| < m}, 

W'"’°°([0,r] X 5) = {fit,y,z) e F°°([0,r] x 5)1 Z“/ e L“([0,r] x 5), |a| < m}, 
and the following notations 


l^m.oo = sup ||Z“/(T)||icx=( 5 ), m > 1, ||/||oo,t= sup ||/(T)||i.x,( 5 ). 

‘ 0<r<t I ^ 0<r<t 

— — |Q;|<m — — 


2.1 General Inequalities 

The following Gagliardo-Nirenberg-Morser type inequality will be used repeatedly whose proof can be 
found in |19) . 

Lemma 2.1 The following products and commutator estimates hold: 

(1) For u,v G L°°([0,T] x 5) fl W^([0,r] x S) with k G N be an integer. It holds for any 0 <t <T that 

f \\{2^uZ''v){T)fdT<\\u\\l^tf \\v{T)\\l^^dT+\\v\\l^t f ||u(r)||^,dT, \l3\ + \iy\=k (2.1) 
Jo Jo Jo 

(2) For I < \a\ < k, 0 < t < T, g G W'=-i([0,r] x 5) nL'^([0,r] x S), f G W'=([0,T] x S) such that 
Zf G L°°([0,T] xS), we have 

f\\[Z-J]grdr<\\g\\l^, r||Z/ll«-irfr + ||Z/||L,t T ( 2 . 2 ) 

^0 -'0 -'0 

(3) For |a| = fc > 2, 0 < t < T, and f,gG >V^“^([0,T] x S) fl >V^’°°([0, T] x 5). Define the symmetric 
commutator as [Z, /, g\ = Z°‘[fg) — Z° fg — fZ°‘g. Then, it holds for any 0 < t < T that 

f\\[Z-,f,g]fdr<\\Zf\\l^, f\\Zg\\l,,.,dr+\\Zg\\l^, f \\Zf\\l,,^.dr. (2.3) 

Jo Jo Jo 

We also need the following anisotropic Sobolev embedding and trace estimates in S whose proof can 
be found in [32]: 

Lemma 2.2 Let mi > 0, m 2 > 0 be integers, f G fl and V f G 

1) The following the anisotropic Sobolev embedding holds: 

ll/lli=.<c'(||v/iu, + ||/iuj||/iu,. (2.4) 

provided mi + m 2 > 3. 

2) The following trace estimate holds: 

< cdlv/IU. + ll/IUJIl/IU,. (2.5) 


with mi + m 2 > 2s > 0. 

The following classical Sobolev inequalities and commutator estimates in hold: 
Lemma 2.3 For s G R., s > 0, it holds that 

< C's(|/|L~(R2)|g|_H-<»(R2) + |5U“(R2)|/|ff«(R2)), 

< C's(|V/|ioo(R2)|5|^a(R2) + |Vg|i=o(R2)|/|^„(R2)), 

|uu|l < |u|i,ook|i, 

where is a Fourier multiplier of {1 + |^P)2. 
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2.2 Estimates involving (p 

We shall derive some estimates of 77 chosen by (I1.19|) and (11.201) in terms of h and v. It shows that rj has 
^ higher order standard Sobolev regularity in S than h. Note that the different choice of cut-off function 
K in (11.201) from the one in |32] allows one to control the norm of rj itself. Assume that A is chosen 
such that dzifioiy, z) > 1 at the initial time and 

dz‘p{,t,y,z) > Co, Vte[ 0 ,r^], ( 2 . 6 ) 

for some cq > 0 and small time T®. With this a priori assumption, one obtains that 

Lemma 2.4 It holds forrj, defined in (I1.19l) . (ll.20l) . that: 

VfceN, 

'ik,lGN, \\d[r,\\H^^s)<Ck\d[h\,_i, 
and moreover, we also have the L°° estimates 


Vfc e N, < C'fc|/i|fc,oo, 

Vfc,/eN, \\d[y\\wK^<Ck\d\h\k,^. 

Proof. It follows from the expressions (11.191) and (11.201) that 

\\yW^=f dzf m,f,z)\H^<[ \h{t,0?dff 

' {0-^\h{t,0\^df = \h\.. 


Hziordz 


< 


(2.7) 

The proof of remaining parts is similar to |31j . so we omit the details here. □ 

Lemma l2.1l and [2.4l vield immediately the following estimates which will be used repeatly, whose proof 
can be found in EH. 


Lemma 2.5 For any m G N, it holds that 


f \\^\\nr.dr < A(l, ll/IU.t + r WfWnr. + \Z^h\ldT, 

JO ^0 Jo ^ 

and moreover, for standard Sobolev norms, we also have that, for 0 < k < m, 

pt p -I pt ^ 

/ ||a^-(7^)llff™-rfT<A(-,||/|U,t + |/rU)-) / + 

Jo Ozif Co Jo ,__r, ^ 


( 2 . 8 ) 


(2.9) 


k=Q 


As a consequence of Lemma 12.51 one has 

Corollary 2.6 For any m > 1 and any sufficiently smooth function f, it holds that 

f WV^ffur^dT < K{-, |/r| 7,0.) f ||V/|||,™dr 

Jo Co * Jo 

+ H-, \\dzf\Ut + |/i|w=.=o) f + \Z^h\ldT, 

Co Jo ^ 

f W^fWlr^dr < A(l, |h| .,00) r ||V^/||?,„dr 
Jo Co * Jo 

+ A(l, || 5 ,/||oo.* + |/i|w=,=») r WdzfVn^-^ + \Z^h\ldT, 
Co Jo ^ 


( 2 . 10 ) 


( 2 . 11 ) 


f \\^‘^V^ffnmdr < A{-,\h\^.^){ /‘||VVll?^™dr+ T ||V/||?,„dr 
Jo Co * l. Jo Jo 


IVVIIL .7 r \Z^h\ldr + IIV/IIL., r \^yZ^h\ldr}. ( 2 . 12 ) 

Jo Jo 2 ^ 


and 
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Proof. These estimates follow from (12.21) . (12.81) and the identity, for i = 1,2,3, 

Z“(9r/) = - [Z“, ^]dj, 

Oz 9^ Oz ^ 

where ditp should be replaced by —1 when i = 3. □ 

Corollary 2.7 For 1 < |q;| < /c, 0 < t < T, and any sufficiently smooth function /, it holds that 

f \\[dt 2“]/f rfr < A(l, ||V/||oo.t + |/i|w?-) r l|V/||^.-i + \Z^h\ldT. (2.13) 

do Co do ^ 

Proof. This corollary follows from (12.2L (12.811 and the identity, for i = 1, 2, 3, 


[df,z°‘]f = [z“ 


Ozif Ozip 


where dtip should be replaced by —1 when i = 3. 


□ 


The following lemma gives the regularity of free surfaces once we gain the corresponding one of the 
velocity. 

Lemma 2.8 For any m G N, e G (0,1], it holds that 


s\Z^h\ldT < elZ^hoW 
2 2 



II + A(|Vy/i|oo,t + ||?^|lw)'“ 



e\Z^h\l)dT. (2.14) 
2 


Proof. The proof of this lemma is similar to Proposition 3.4 in [32] except for the time derivative 
involving. By taking time derivative in the boundary condition (11.321) and following the same argument 
in m, we can prove this lemma. The details are omitted here for brevity. □ 


Since the Jacobian of the change of variable (11.2411 is dzip, it is natural to use on S the following 
weighted scalar products: 


fgdVt, with dVt = dzip{t,y, z)dydz. 


(2.15) 


With this notation, one has the following integration by parts identities for the operators df : 
Lemma 2.9 For any sufficient smooth function f and g, it holds that 


where N = {—dih, 


[ dffgdVt = - [ fdfgdVt + [ fg-N.dy, 

J S Js Jz—0 

[ dffgdVt = dt ( fgdVt- f fdfgdVt- f 

J S S ^ S 'J z - 

-d2h,iy. 


* = 1,2,3, 

(2.16) 

fgdthdy, 

(2.17) 


0 


The following Korn’s inequality is also needed(see [32] for the proof). 

Lemma 2.10 Assume that dzip > cq and ||V:p||loo + ||^ for some cq > 0, then there exists 
Aq := A(^), such that for every v G F[^[S), it holds that 


and 


||Vz;|| 


||V^;||i.(5) < Ao / IV^^vl^dVt, 

Js 

i2(5) < Ao^^ \S'^v\^dVt + ||'c|li2(5)^, 


(2.18) 

(2.19) 


where S^v = ^{V^v + 
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3 Equations for higher order conormal derivative terms 

In this section, we derive the equations for {Z°‘v, Z°^p, Z°‘rf) involving high order conormal derivatives. 
As explained in [32] , the commutators between Z“ and df + u • or in the equations yields a loss 
of i derivative due to the lower regularity of ip when cr is not a given constant independent of e. We shall 
introduce the Alinhac good unknowns [I] to yields some nonlinear cancellation in the equations so that 
the commutators are all the lower order terms. 


3.1 A commutator estimate 

In order to perform higher order conormal estimates, we first derive the equations satisfied by {Z°‘v, Z°‘p). 
We thus need to commute the conoraml vector field with each term in the equation (ll.llll) . The 
following commutation relations are needed, which has been shown in [32] except for Zq = dt- 
It is easy to check that for * = 0,1, 2, 3, and a smooth function /, one has 


= dfZ^f - dUdfZ^rj + Cnf), 

where the commutator Cf{f) is given for a 0 and i 3 by 

c“(/) = c“i(/)+q2(/)+c“3(/), 


(3.1) 


(3.2) 


with 


c“i(/) = -[z“,|^,a,/], 

Oz^ 

c?,u} = 


For i = 3, we only need to replace dip by —1 in the above expression. The following commutator estimates 
hold: 


Lemma 3.1 For 1 < |a| < m, i = 0,1,2,3, one has that 

f l|C“(/)f dr < A(l, |/i| ..oo + llV/IUcoo) f ||V/||^_, + \Z^-^h\ldT. (3.3) 

Jo Co * Jo ^ 

Proof. The proof is similar to m which involves only spacial norms, since the role of Zq = dt is the 
same to tangential derivatives. The details are omitted here. □ 


3.2 Equations satisfied by Z'^p, Z'^'q) 

We now derive the system satisfied by Z'^v, Z°^p, Z°‘t]. Alihnac’s good unknown are defined as 

V“ = Z°‘v - and - dyz°‘r]. (3.4) 

It should be noted that the following computation work only in a neighborhood of the boundary, but, for 
notational convenience, we write it in a general form without confusion. 

Lemma 3.2 For 1 < |a| < m, V° and Q°‘ satisfy 

div^V°^ + —dtQ°‘ + —u • 

IP IP 

= Z‘^rj{dt{—)dtp + dti — ) ■ VV) - C“(d) - —Q{p) - — ■ C^{p) 

7P 7p 7p 7p 

=--K, 


(3.5) 
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and 


Q{dt + V ■ V^)F“ + V‘^Q“ - 2 nediv‘^S^V‘^ - \eV^div'^V°‘ 

= Z°‘r]d‘^Q{dtv + V ■ V^z;) + gZ°‘7]{d‘^v ■ V‘^)v - [Z“, g](dfv + v ■ V^v) - C°^{p) 
- gC°‘{T) + 2nediv'^£°‘{v) + 2g,eV°‘{S^v) + AeV^C“(d) + \eC°‘{div^v) 


=: 7^' 


Ml 


where the commutators Cq {p), C°‘{p), C°^[d), C°‘{T) and£°^(v) admit the estimates: 

f ||(Co“(p),C“(p))f dr < A(l, \h\^.^ + llVpIl .,=») r ||Vp||?,._. + \Z^-^h\ldr, 

Jo Co ^ * Jo ^ 

f ||C“(d)f + ||£:“(z;)f dr < A(l, \h\^.^ + |lVz;|| ., 0 .) f ||Vz;||?,™-. + |Z—^dUdr, 
Jo Co ‘ * Jo ^ 

f ||C“(r)f dr < A(l, IKz;, Vz;)|| i,o„ + |d|^..») f \\v\\l^ + \\Vv\\\^., + |Z™d|i.dr, 
Jo Co * ‘ Jo 


(3.6) 

(3.7) 

(3.8) 

(3.9) 


and 


II (dzz;‘" £:“(?;), 27“ (5‘^z;), V“C“ (d), C“ (dzz;‘"z;)) f dr 


(3.10) 


< A(l, |d|^3,=„ + II Vz;|Ui,=. + f £^11 V\||^._, + II Vz;||^_i + £^1 v,Z™-id|idr. 

Co * * Jo ^ 

Proof. It follows from (13.11) that 


Z“div‘^?; = div‘^Z“z; - • V'^Z“r; + C“(d), 

Z“(^5rp) = ^{dtZ^p-dtpdtZ-p + C^p)) + [Z^,^]dfp, 

2 :“(— . V» = — • (V’^Z^^p - dtpV^Z°‘p + C°^{p)) + [Z“, —] • V'^p, 

7P 7p 7P 

where C“(d) = X)i=i^f(^j)i C°‘{jp) = (C“(p),C2 (p),C3 (p))‘. It follows from Lemma (XU that 

f ||C“(d)f dr < A(l, \h\^.^ + llVzzIUj.oo) f ||Vz;||^_, + |Z™-id|idr, 

Jo Co * * Jo = 


||(Co“(p),C“(p))f dr < A(-, \h\^.^ + llVpIl 

Co * 


W! 


llVpIl 


7 m —1 L |2 


d|ldr, 


(3.11) 

(3.12) 

(3.13) 


which proves (13.71) and a part of (13.81) . Adding (I3.11I) - (I3.13I) and using (ll.31|) ^. one obtains (13.51) . 
Note that 


Z“(£<(5fn + V ■ V^zi)) = gZ°‘{dfv + v ■ V'^n) + [Z“, g]{dfv + v ■ V'^v). 
Rewrite the transport operator df + v- as 


(3.14) 


df = dt + VyVy + V,d,, (3.15) 

where 

14 = ■ N - dtp) = -^{v ■ N - dtp), (3.16) 

dzP dzp dzp 

and N(t, y, z) defined as 

N(t, y, z) = (-5i77(t, y, z), -d 2 p(t, y, z), 1)‘. 

is an extension of the outward normal vector to the boundary. Henceforth, we will use it and the outward 
normal vector to the boundary without confusion. 

It follows that 


Z°‘{df + V ■ V‘^)u = (df + V ■ V^)Z“?; - d^v{df + v ■ V^)Z°‘p + C“(T) 


(3.17) 
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where the commutator C°‘{T) is given by 


6 


c“(r) = 


la 

i 5 


and 



Finally, for the viscous term, it follows from (EU that 

Z°‘div^S'^v = dW°‘Z°‘S'^v - ■ V^Z°‘r] + V°‘{S'^v), 

where 'D°‘{S‘^v)i = C‘^{S'^v)ij. Note that 

= S'^Z°‘v - - V‘^Z“r; (g) 9^ + E°‘{v), 


(3.18) 


(3.19) 


with {£°‘v)ij = Cf{vj) +C“(?;i). Therefore, it follows from (13.181) and (13.1911 that 

Z“div‘^5'‘^z; =div^S'^Z“u - div‘^(a> g) + V^Z“7? ® 

+ div‘^£:“(z;) - + 2iieV°{S'^v). (3.20) 

Similarly, one has 

Z“V^div^u = V^Z“div^z; - 9^div‘^uV^Z“7? + C“(div‘^u) (3.21) 

= V^div^Z^u - V^(a> • V^Z^r;) - a^div^uV^Z“77 + V“C“(d) + C“(div‘^u), 

where (|3.11l) has been used in the last equality. Collecting (I3.14|) . (13.1711 . (I3.20F (13.2111 and using (11.3111 ^ 
lead to (13.6F Then Lemma [3.11 yields that 



2 


and 




Co ‘ ‘ 


+ ||Vu||^^-i +e2|VyZ™-i/i||dr. 


It remains to estimate the commutator C°‘{T). One first needs some estimates of Vz and 14 which will 
be used later. It follows from the definition of Vz and 14 in (13.1611 that 



(3.22) 


(3.23) 


Using m and Lemma [211 one has that 



(3.24) 






















20 


It follows from Lemma [23] that 


[ WZV^Wl^m-.dT < [ + 

Jo Jo 'JzV 

1 /"* 

< A( , IjrJ^llyyl.OO + |ft,|yy2,Oo) / \\{y ^ , Z V z)\\-^'m, -2 + \ 
Co * Jo 

1 /'* 

— A( , ||l'||yyl,oo + |/l|yy2,oo) / ||r>||^m-l + \h\-j^mdT. 

Co * Jo 


?m-2j,|2 , II 2 


'/ill dr+11 


i-2dr 


As a consequence of Lemma [1+] and the above estimates for Wz, 14, one gets that 

r ||C“(r)f dr < A(l, ||(u, Vz;)|Uc=o + \h\^ 2 ^) f \\v\\l,^ + || + \h\l,^dT. 

Jo Co ‘ * do 

Therefore, the proof of Lemma [321 is completed. 


(3.25) 


□ 


3.3 Boundary conditions for {Z^^v, Z^^p, Z^r]) 

In the following, we derive the boundary conditions for {Z°^v, Z°^p, Z°‘r]). Note that the only interesting 
case occurs when = 0, since Z'^v = Z°‘ri = Z°‘p = 0 with 03 0 on the boundary. We start with the 

dynamic boundary condition. 

Lemma 3.3 Let 1 < |a| < to with 03 = 0, on the boundary {z = 0}, it holds that 

{2peS^V°^ + Xediv^V^)!^ + Z°‘h{2ped‘^S’^v + Xed^div‘^v)N - aVy ■ Z“ ( ) N 

V Vl + 

- Z>N + ( 2peS^v + Xediv^v - aVy ■ ^ -{p- Pe) ] 2:“N = C“(B), (3.26) 

V Vl+|Vy/lp / 

where the commutator C°‘[B) satisfies 

f |C“(i?)|i2dr < A(l, \h\y,2^ + llVulUj, 00)^2 ^ |V,/i|l,„_,dr. (3.27) 

do Co * * do 

Here and henceforth, ()^ denotes the value of function on the boundary z = 0. 

Remark 3.4 It should be emphasized that the (13.261) holds only locally on the boundary. But, this is 
enough since one can assume that the solution {g,v) is supported compactly on a neighborhood of the 
boundary. While, for the interior domain the boundary condition will not be involved in the conormal 
estimates. 


Proof. Applying the operator Z“ to (11.331) and using (I3.11|) . p.l9l) . and the Alinhac good unknowns, 
one obtains (|3.26l) with the commutator C°‘{B) of the form 


C“(R) = [Z^,p-pe + cTVy 


V„/i 


" v/l + |V,/iP 

In order to control C°‘{B), we first notice, from (I1.33L that 


— 2pL£S'^v — Aediv‘^v,N] — 2fj,££°‘{v) — AeC“(d). 


p- Pe + crVy 


Vyh 


yr+rw 


= 2ps{S'^vn) ■ n + Xediv‘^v. 


which implies that 


/ |C“(R)|i2dr<£' / |[Z“,n(5^u),N]|i.(R.) + |£:“(«)| 
^0 ^0 


L2(R2) 


|C“(d)| 


i2(R2) 


dr. 


(3.28) 
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Here and henceforth, H := Id — n 0 n denotes the tangential vector field. It follows from Lemma 1271] that 

Jo Jo 

< A(l, \hU^ + ||Vi;|Ui,^) r + \Vv%m-^dT. (3.29) 
Co * Jo 

To estimate £°'(v) and C“((i) on the boundary, one needs only to control f* |C“(uj)|| 2 (R 2 )d'r. Note that 
Cf^{vj) = 0 since we only consider the case that = 0. Then, using the similar argument as in the 
proof of (EJl), one gets that 

f \Cnvj)\^dT = A(l, + llVull i.oo) f + |V,h|^^_,dr. (3.30) 

Jo Co ‘ ‘ Jo 

Substituting (13.2911 . (13.3011 into (I3.28|l . we prove (13.2711 . Therefore, the proof of Lemma 13.31 is completed. 

□ 


Lemma 3.5 For any 1 < |a| < m such that ao = 0, on the boundary {z = 0}, it holds that 

dtZ^^h - ■ Z“N - • N = C“(h), 

where the commutator [h) satisfies the estimate 


\C°{h)\l2dT < A( —, ||u||wi.=o + ||Vu||oo.t + \h\y^2^) [ 
Co * * Jo 


+ V'. 


b\2 

ylW” 


Proof. Applying the operator to the kinematic boundary condition (11.3211 ^ shows that 
dtZ^h - ■ Z“N - • N = -[Z“, v^, Vyh] + := C“(/i). 

Uz^ 

It follows from Lemma [O that 

f \C^{h)\l2dr < A(l, Hull 1,=. + llVulloo.t + t \h\l,^ + 

Jo Co ‘ Jo 

which complete the proof of this lemma. 


(3.31) 


(3.32) 


(3.33) 


□ 


4 A priori Estimates 

The aim of this section is to derive the a priori estimates in Theorem 14.11 which is a crucial step to prove 
Theorem 11.51 We drop the superscript e throughout this section for notation convenience. 

Theorem 4.1 (A priori estimates) Let m be an integer satisfying m > 6, e G (0,1], cr G [0,1], and 
{p,u,F) be sufficiently smooth solution, defined on [0,T®], to p.lll . (11.511 . (11.611 . (11.811 . (11.911 with the 
initial data satisfying (11.3811 and (11.4611 . Then it holds that 

l£'o(y,^)|exp(- f \\div^v{T)\\L--dT) < g{y,z,t) < \go{y,z)\exp{ f ||dzt;'^u(T)||LoodT),Vt £ [0,T""]. 

Jo Jo 

(4.0.1) 

Furthermore, there exists a time Ta > 0 independent of a and e such that the following a priori estimate 
holds: 

0m{t)= sup {l + \\{p,v,h){T)\\x^)+ [ ||Vp(t)||^„_i + ||Ap(r )||^2 + ||Vu(r)||^„_idr 
0 <r<t Jo 

+ e f ||Vu(r)||^™ + ||V\(T)||^™_ 2 dT + e2 /" ||V\(T)|l^,„_idT 
^0 -'0 
<A(C'o,Co), vte [0,mm(r„,r,)] 

where Ta depends only on Co, do and Co- 


(4.0.2) 
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The main difficulties in the proof of Theorem 14 .1 1 are to derive higher order conormal energy estimates 
for (p,v,h). To overcome these difficulties, we will apply the Alinhac good unknowns 

= Z°‘v - Q“ = Z> - a ^ 0, (4.0.3) 

and for a = 0, we denote = v and = p. As emphasized before, the good unknowns are chosen so 
that the corresponding commutators in the equations for 14“ and (5“, which are computed in the previous 
section, are all lower order terms. The control of these good unknowns and Z'^h, which can be obtained 
by the standard energy method, yields a control of v, p and h. Set 

\\v^[t)r= E = E E 

|a|<m |Q;|<m |Q;|<m 


Throughout this section, we assume the following a priori assumptions hold: 


0 < — < g{t) < 4Co, > Co, 
40o 




1 

Co’ 


Co 


-5^pU=o> Vte [0,T^]. (4.0.4) 


Thus, by using the assumption (|4.0.4|) . one has that 


< \\v^{t)r + E \\dtv{t)z-vit)r < \\v^{t)r +a( 1 , \\vv\M\h\i^, (4.0.5) 

|a|<m 

umiir. < \\Q^{t)r+ E \\dtp{t)z-p{tw < ( 4 . 0 . 6 ) 

|a|<m ° 


Remark 4.2 As explained above, throughout this section, it is always assumed that {g,v,h) has compact 
support in or in a vicinity of dS, and the system (11.311) holds in or S because one can realize 
this assumption by multiplying a cut-off function and at the cost of some source termsfwhich are easy to 
estimate in the conormal Sobolev space because at least one derivative is applied to the cut-off function) 
appearing in the right hand side of (11.311) . 

Throughout this section, we shall work on the interval of time [0,T®] such that (14.0.41) holds. And 
we point out that the generic constant C may depend on ^ x jn this section. Since the proof of 
Theorem 14.II is very complicated, we divide the proof into the following subsections. 


4.1 Basic energy estimate 

First, we give the basic energy estimate which corresponds to the physical energy. Since it is not necessary 
to use the local coordinates to prove the basic energy estimate for we use the Euler coordinates 

here. Precisely, the following lemma holds: 

Lemma 4.3 For a smooth solution to uni), (usi) and (HU, it holds that 

f dx-\-pe\ilt\ + + [ [ 2p,£\Su\'^-\-Xe\divu\^dxdT 

dOt 2 7-1 Jo Jn^ 

< [ EoMq + +pg|fJo| + (t|So|. (4.1.1) 

J Oq 2 T 1 

where |flt|, |St| denote the volume of Fit o,nd surface area ofEt, respectively. 

Proof. Multiplying (ll.ll) „ by u and integrating over Fit, one gets from integration by parts and the 
kinetic boundary condition (HU that 




pdbfudx-\-2pLe / \Su\'^dx -\- Xe / |di 

I; ^ Clt 'J r2t 


^dx 


+ 



2pLeSun — Aedivun). udS. 


(4.1.2) 
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where we have used the mass equation (ini)i and the transport formula. The dynamic boundary condition 
m yields that 



2^eSun — Aedivun) • udS 



aHn) ■ udS 


Pedivudx — a / Hn-udS. 


(4.1.3) 


Using the mass equation (ini)i and the transport formula, one obtains that 


— / pddvudx + / Pedrvudx 

J r2t r2f; 


d 

dt 





(4.1.4) 


Note that the mean curvature of the free surface is give by (inD. Let Et be determined by a local 
coordinate system x = a;(si, S 2 , t)^ {si, S 2 ) € U so that its Riemannian metric given by 

ds = gapdx°‘dx^, 


where U is an open set, gap = Xa ■ xp with Xa = dsc^x and the convention summation over the repeated 
indices has been used. It is well known that the Laplace-Beltrami operator can be written as 

Astit) = ^^s^^{^/gg°‘^^,|3), a,/3=l,2, (4.1.5) 

V9 

with g = detgap and (g“^) = {gap)~^■ Then, one gets that 


— cr / Hn-udS = —a f g ^ ds<^{-\/gg°‘^d^px) ■ udS 
Jst Jut 

= -cr ds<^{y/gg°‘^dsiix) ■ xtdsids2='^ / y/gg""^dtgapdsids2 
Ju ^ Ju 

= [ y/gdsids2 = CF^ [ dS. 

dt Ju dt Jqq 


(4.1.6) 


Substituting (I4.1.3l) - (j4.1.6p into (14.1.21) . one gets that 


-^( / -pu^ + dx + Pe\^t\ + alAttl) + f 2/xe|S'up + AeldivuPda; = 0. 
dt 2 7-1 

Integrating the above equation, we proved (14.1.ID . Therefore, the proof of this lemma is completed. □ 


4.2 Higher order conormal energy estimates 

Next, we perform the standard energy estimates on the Alinhac good unknowns U“ and Q" solved by 
( 1 ^ . dm). However, a new difficulty shall be overcame when the effect of surface tension((T ^ 0) is taken 
into consideration. It seems difficult to bound the boundary term involving the highest order derivatives of 

yh\ 

condition (I3.26P and the equation (13.6D to turn this boundary term back into the volume integral. On the 
other hand, since (I1.62D is invalid for compressible flows, || V^u||^m_i(ir will appear in the estimate 
of U“ and ( 5 “, when one directly apply the trace estimates to terms involving Jq |(9zu)^|^m_i(ir(cf. 
Lemma EH- Fortunately, with a there, this term can be bounded by the dissipation terms in the 
estimates for normal derivatives and div'^u in the next steps. Precisely, one has the following lemma. 


crVy 


y/l+Jv 


due to the less regularity of h, ^Acf. p.BdP b We will use the third equation in the boundary 


Lemma 4.4 For any t £ [0, T^] and m > 5, it holds that 

||U™(t)f + ||Q-(t)f + \hit)\l,^ + + e f IIdr 

do 

< A(l, ||U™(0)f + ||Q’"(0)f + ||u(0)||?,™ + alVyhml,^ + |M0)|?i^) + Aod f || Vp||?,_.dr 
Co Jo 

+ AoSe^ [ \\V\\\l^m-idT + CsAaoit) [ A(L;„(t)) + ||Vu||^.„_i + e|Z™/i|idr, (4.2.1) 

Jo Jo ^ 
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where Aq = A(^,C'q), A^oit) = A(^,C'o, Q(t)) anrf 

Q{t) := Vz;)||y^;i,.» + ||(p, Vp)||y^;i.=o + \h\^3,^ + ei\\d^^v\\oo,t, (4.2.2) 

y„(i) := ||(F-,Q-)(t)f + ||V(p,^)(i)||?,_. + \{h,y^Vyhm\l^ + WA^pmh 

+ e\\\7{p,v){t)\\l^^_, + e\\A‘^p{t)\\l^ 2 . (4.2.3) 

Proof. The case for m = 0 is proved already in Lemma 021 Assume that (14.2.111 is proved for fc < m — 1 . 
We shall prove that it holds for fc = m > 1. Multiplying the equation dSH) by and integrating over 
the domain 5 yield 


_d /■ 1 

dt Jc 2 


p|y“|dVt + J 2ps\S‘^V°‘\^ + Xe\dw^V°‘fdVt 

(2/re5'‘^P“N +AedTO‘^P“N-Q“N)P“d2/+ [ Q°^div'^V‘^dVt + [ n‘^V°^dVt 

Js Js 


'2 = 0 
=:B + J+ [ 


(4.2.4) 


where we have used the integration by parts and the boundary condition (ll.32|) . 
For the boundary term B in (14.2.41) . it follows from (13.261) that 


Bdr = (T 




Vyh 


0 .^2 = 0 


■x/TTrW 


)Bi-V°‘dydT+ f [ d‘^pZ°‘hN ■V°‘dydT 
Jo Jz=0 


f [ ( 2pLeS^v + Xediv'^v — aVy 

Jo Jz=0 \ 


Vyh 


Vl + |V,h|- 


- (p - pe) Z°‘N ■ V^dydT 


t t A 

- f f Z°‘h{2ped‘^S‘^v + Xedyiv'^v)NV°‘dydT + [ [ C°‘{B)V°^dydT =:^Ki, (4.2.5) 

./0 J z—0 Jo J z—0 

which will be dealt term by term. First, we deal with the most difficult one Ki. The kinetic boundary 
condition (13.311) yields 


A'l = a 


VyZ^ 


){dtZ°‘h-v^ ■Z‘^N-C°‘{h))dydT. 


0 .^ 2=0 


Note that 


A/yh 


■,/iTWJiW 

Z'^Vyh Vyh{Vyh,Z°‘Vyh) 


(4.2.6) 


./iTWJW ^l+\Vyh\^ (l+|V,h|2)i 


-C“(5), 


where C“(5) = [ 2 :“,Vyh, 


]. It follows from the integration by parts that 


‘ '• Z“Vy/l Vyh{Vyh,Z°‘Vyh) 


0 Jz=0 "^l+WJW {l + \yyh\^)i 


)dtZ°‘hdydT 


^ _ VykiVyh^Z^^Yh ■ da-Vyhdydr 


0 Jz=oWT+W?^ (l + |Vj;h|2)f 


-f-l (- 

Jo dr J^ 


Z^Vyh\^ 


— o 


< —a 


=0 2v^l + |Vy/i|2 2(1 +|Vj,/i|2)l 

* '■ , \Z‘^Vyh\'^{Vyh,dtVyh) , i\{Vyh,Z‘^Vyh)\'^{Vyh,dtVyh)^, , 

(-5--5- }aym 

10 Jz=o 2(l + |V,,h|2)l 2(l + |Vj,h|2)l 

'•* '• {Vyh,Z°‘Vyh){dtVyh,Z°‘Vyh) 


0 Jz^o 


(l + |V,h|2)l 


|Z“V,h(t )|2 


-dy + cr ( 
Jz^O 


dydr 

Z^Vyhol^ |(Vyho,Z“Vj,ho)|2 


=0 2(1 +|Vj,/i(f)|2)f A=o 2v'l + |Vyho|2 2(1+|Vj,ho|2)§ 


)dy 


+ Aoo(t)cr [ \Vyh\l^mdT, 
Jo 


(4.2.7) 
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and 


^ _ 

0 ^.=0 Vl+ ~ (l + |V,,/iP)i 


Vy/l(Vy/l, Z°'Vyh) ^ 


= CT 


< CT 


(- 


Vyh{Vyh, Z°‘Vyh) 

0 Vz^oVl + l^y^P (l + |Vj,/i|2)i 

Vy/l(Vy/l, Z°‘Vyh) 


)-Vy(i;^-Z“N)dydT 


Z°‘\7yh 

(- 


+ CT 


0 A=oVrTrv7^ (i + iVy/i|2)i 
* '■ Z“Vy/l Vyh{Vyh,Z°‘Vyh), 


^''')-ivt-^y)Z‘^^yhdydT 


, ( , _ 

0 7^=0Vl + |Vy^P (l + |Vy/i|2)5 


(Vy?;'’ • Z“N)dj/dT 


< cr 


< • V,(|Z“V,/i|2) . V,(|(V,/i,Z“V,/i)|2) 


0 Jz = 0 2a/ 1 + |Vy/l|^ 


2 (l + |V,/i|2)^ 


-dydr 


+ CT 


+ CT 


/o 

rt 




=0 (1 + |V,AP)4 

Z°‘V,,h 


0 Jz=o\/T+\%W (l + |V,,/i|2)f 


Vyh{VyKZ-Vyh) ^ ^ . _ Z“V,/.)dydr 


<Aoo( 0 ^ [ |Vy/i|^mdr. 
Using Lemma 12.11 one has that 


and 



(4.2.8) 


(4.2.9) 

) [ \Vyh\l^mdT. 

JO 

(4.2.10) 


Therefore, from (14.2.101) . (I3.32|) and trace estimate (12.5L one obtains that 

a / / VyC°{S){-v'’■ Z°N-C°{h))dydT 

Jo Jz =0 

< aA^it) |V, •C“(5)|i.|' + |C“(h)|2dr 

< Aoo(t) [ \h\l^m + ||t;||^™-i + a\Wyh\j^m + ||Vr;||^™-idr 

Jo 

For /p Vy • C°‘{S)dtZ°'hdydT, if Z°‘ ^ 9™, then (14.2.10|) yields immediately that 

cr [ f VyC’^{S)dtZ°‘h <A^{t)a f |Vy/i|^™dr. 

Jo Jz =0 Jo 

For Z“ = 9™, it follows from integration by parts and (14.2.101) that 


(4.2.11) 


(4.2.12) 


f f VyC°‘{S)d^+^hdydT = a [ f C°‘{S)dJ^+^VyhdydT 
Jo Jz—0 Jo Jz—0 

< [ (7C°‘{S)dTVyhdy + a f C°‘{So)dTVyhody + a [ [ dtC^{S)dTVyhdydr 

J Z—0 J z—0 Jo J z—0 

< Ao[|Vj;/l(0)|^m_l + cr\yyholl^m] + ■^ y 

+ Aoo(t)cr f iVyhlli^mdT, 

Jo 


\Z^^yKtr 

z=o 2(1 + |Vy/i(f)|2)i 


dy + Ao(t\V yhlh^m-i 


(4.2.13) 
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where one has used that, for m > 5, 

m —1 

\c^{sm\h < E 


Vl + lV./iP 


lL 2 


< 


1 


Vi + IVy/iP 


m —1 


li ,oo 


i9r“v,/.ii.+ E 


[^]+i 

2 


x/T+rw 


\h 


Co "■ 


It should be reminded that the estimate of 


0 Jz^o Wl+\Vyh\^ (l + |Vj,/i|2)i 


(4.2.14) 


(4.2.15) 


is a subtle issue. If one uses directly the integration by parts, the term |, which is difficult 
to estimate, must appear due to the boundary term Z^Vy- ZVyh with |/3| = m— 1 in C°'{h). Therefore, to 
deal with this boundary term, we will take full advantage of the boundary condition (13.261) and equation 
to turn it into a volume integral. In fact, the dynamic boundary condition (13.261) yields that 


-aVy •Z“(- 


^ -=) - d^pZ^h 


" Vi + |v,/iP' 

= Q“ - {{2peS^V°‘ + Aediv‘^y“)N)3 - {Z°‘h{2ped‘^S’^v + Ae5^div^u)N)3 + C“(B)3. (4.2.16) 

where (•)3 denotes the third component of a vector. Thus, it follows from (13.271) and trace estimate (12.51) 
that 

-af f VyJ "l Z^vl ■ ZVyhdydr 

Jo Jz=o \VTT\vJW (i + |Vj,/i|2)i ) ^ 

f [ Q°‘Z^vl-ZWyhdydT- [ f ({2peS‘^V°‘ + Xedw^V°‘)N) Z^vl'ZWyhdydT 
0 Jz =0 Jo Jz =0 ^ 2 3 

+ Aoo(t) \v%^-.dT^ ' \h\l,^ + |C“(i?)3|i. + <j\Vy • C^{S)\l.dT^ 


< 


< 


/o Jz^Q 

rt 


Q^^Z^vl ■ ZVy 


hdydr-J ({2pLeS^V^+ \ediY^V^)'N^^Z^vl-ZVyhdydT 

rt 


S f S^\\'X/‘^v\\‘^m~ldT CsAcait) f + ll'l’ll^m -1 + II Vr>||^m -1 + cr| Vy/lj^rndT, ( 4 . 2 . 17 ) 

Jo Jo 

where |/3| = m — 1. To estimate the boundary term fg ■ ZVyhdydr in (I4.2.17L using 

integration by parts, we note that 


[ f Q°‘Zf^vl ■ ZVyhdydr 

Jo Jz=0 

= [ [ d^Q°‘Z^vy ZVypdVrdT+ j f Q°‘d‘^{Z^Vy ■ ZVyy)dVrdT 

Jo Js Jo Js 

rt 


( 4 . 2 . 18 ) 


/o JS 
rt 


< f [ dtQ‘^ZPvyZVyydVrdT+ [ ||Vu||^,„_idT + Aoo(t) / ||u||^„._i + ||Q™f dr. 

Jo Js Jo Jo 

Using the equation (13.61) . one obtains that 

[ [ d‘^Q°Z<^Vy ■ ZVypdVrdr 

Jo Js 

<- [ ( Q{dt+vyVy + VA){V^)gZ^vyZVyydVrdT ( 4 . 2 . 19 ) 

Jo Js 

+ [ f i2yediv^S‘^V^ + XeV‘^div^V°‘)3Zf^vy ZVyydVrdr + Aooit) f ||7^^||||v||«n.-ldr, 

Jo Js Jo 
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where \(3\ = m — 1 and ||77.^|| will be estimated later. Integration by parts, and using 14 = 0 on 2 ; = 0, 
one gets that 

- [ [ Qidt +VyVy + V,d,){V‘^)3Z^Vy ■ ZVyTldVrdT 

Jo Js 

= - f g{V^)3Z^vyVyijdVt+ [ goiV^)3Z^Vy{0)-VyrjodVo 

Js Js 

+ Aoo(t) ri|b^“(r)f+ |Kr)||?,™dr 
^0 

-hi + Ao[||^“(0)f + \\v{0Wh^-.] + A^{t) + \h\l^dr. (4.2.20) 

For the second term in H4.2.19|l . it follows from the integrating by parts that 

[ f ( 2 /rediv^S'^y“ + XeV'^div‘^V°‘) 3 Z^Vy ■ ZVyj^dVrdr 
Jo Js 

= - [ f Ate(V‘^I/^ + • V‘^{Z^Vy ■ ZVyT]) + Xedw‘^V°‘d‘^iZ<^Vy ■ ZVy'q)dVrdT 

Jo Js 

+ I I (^{ 2 gieSW°^+ Xediv'^V°‘)N'^^Zl^vl-ZVyhdydT 

<Ss [ \\VV"^fdT + CsAao{t) [ e||Vn||^™_x 

Jo Jo 

+ I I ((2MeS'‘^V“ + Aediv‘^fo“)N)^Z^n^-ZVyM 2 /dr. (4.2.21) 

Combining (I4.2.17ll - (j4.2.21ll and noticing the cancellation between the worst boundary terms in (I4.2.17|) 
and (14.2.2111 . it follows from (14.0.511 . for |/3| = m — 1, that 

-u/7 VyJ . ZVyhdydr 

Jo Jz^o {y^T+WJW {l + \yyh\2)i j y 

-hi + + e\\VV^\\lm-^dT 

+ Ao[|ll/™(0)f+ ||u( 0 )||^„._i]+C' 5 Aoo(t) f Y^ + \\Vvf^^-, + \\nU\YmdT. (4.2.22) 

Jo 


For the terms of the form Z^VyZ"'Vyh with |/3| + I 7 I < to, 1 < |/?| < to —2 and the term 




^ .. 

in C°‘{h), it follows from the integration by parts. Lemma [2.11 trace estimate (12.511 and (14.0.511 that 

J« L=a \yTTfvv7 (1 + |V,AP)* y y y J « 

= af [ { ^ (z^^yZ^Vyh + (M 

Jo Jz=o VV 1 + |Vy/i|2 (1 +|Vj,/r|2)i j yy y 

< A^[t) f IIVz;||?,_. + ||y™f + + a\yyh\l,^dr. 

Jo 


Z“/iN 


-Z“/iN dydr 


(4.2.23) 


Substituting (I4.2.7I1 - (I4.2.13I1 . (14.2.2211 and (14. 2. 2311 into (I4.2.6L one obtains, for Z“ ^ 9™, that 


K 3 < Ao(a|V,/io||i^ + ||^’"(0)f + lk(0)||^_x) -a 


|Z“V,/l (<)|2 


-dy 


.=0 4(l + |V,/r(t)|2)f 
+ ^l e|fo“(<)|"dVt + C'7^'l|V"^;||?,™-i+e||Vl/™||?,_,dr 

+ CsA^it) [ Ym + ||Vn||^™_i + ||7^^||r™dr, (4.2.24) 

^0 
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and for Z“ = 9™, it holds that 


K, < Ao[|V,h(0)|^_, + a|V,ho|^™ + ||F™( 0 )r + lk( 0 )||^™-i] + Ao|V,h|^_, 

/o 


— a 


-dy + K^S [ e^VMn^-^ + 4^V^\\^^^^dT 
Jo 


+ 


=0 4(1 + |V,M 0 P)^ 

^ J g\V‘^{t)\^dVt + CsAooit) J Y^ + \\Vv\\\m-^ + \\nit\\Ymdr, 


where the remaining term ||7^^|| in (14. 2. 241) and (I4.2.25|) will be estimated later. 
It follows from (I3.31L (13.321) and (12.51) that 


K. = 


f [ d^pZ°‘h{dtZ°‘h-v'^ ■Z°‘N-C°‘{h))dydT 
Jo Jz=0 

If dtp{t)\Z^Ht)\^dy-l f 9^p(0)|Z“M0)Prfy 


fz^O 


fz^O 


<-k 


f [ dtdtp\Z°‘h\^dydT - [ [ d’^pZ°^h{v^ ■ Z^^N + C°‘{h))dydT 

Jo Jz^o Jq Jz^O 

d'^p{t)\Z°h{t)\'^dy-^ [ d‘^p{0)\Z°‘h{0)\‘^dy + Aoo{t) [ |/i|^^dr 

:=0 ^ Jz^O Jo 

f [ dtpv^ ■ + dtpZ‘^hC‘^ih)dydT 

Jo Jz^o ^ 


/O Jz^O 

<\l dtp{t)\Z<^h{t)\^dy-\ f 9rp(0)|Z“M0)Prfy 

^ J z^O ^ Jz^o 

+ Aoo(t) f ||Vr;||^m_l + + Ilfll^m-ldT. 

Jo 

Since the boundary condition (11.331) yields immediately that 

p — Pe = {2peS‘^v + Xediv‘^v)n ■ n — crVy 


A/yh 


x/T+rw 


This, together with Lemma |2.31 and trace theorem, implies that 


i^3<Aoo(t)e / 


t \ 2 / /** 

Q:TVT|2 i I / \/'i/'a\b\2 


\iv4 \idT 


<5e [ \\VV'^fdT + CsA^{t) [ e\Z"^h\ldT. 

Jo Jo " 

Direct calculations yield immediately that 

Ki< A{—,\\v\\w2,^^t + \\Vv\\wi.^^t)e^\\d^zv\\^^t£^ [ |/i|«m|y“|dT 
Co Jo 

<6e [ WA/V'^fdr+ CsAo.it) [ i\h\l^m + \\V^f)dT. 

Jo Jo 

It follows from (|3.27|) that 

K, < £ |C“(i3)U.|^“lL=rfT < ' (^J* ||r“||||VI/“||^ ' 

< Ao.it){e^Jl ^+ \Vyh\l^.,dry[j'^ l|Vy“||||r“|| + ||r“f}" 

<5e f ||VI/'"fdr + fc2 [ \\v4\\^,n-idT + CsAo.it) f Y^ + \\Vv\\^m-idT. 
Jo Jo Jo 


(4.2.25) 


(4.2.26) 


(4.2.27) 


(4.2.28) 


(4.2.29) 
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Therefore, from (I4.2.24I1 - (I4.2.29L one obtains the boundary estimate for Z°‘ ^ 9™ that 


f Bdr <Ao[|/io|?i™ + a\Vyh^\lr. + ||F™(0)f + ||z;(0)|| 
Jo 




dtp{t)\Z‘^h{t)\^dy 


fz^O 


|V,Z“/r(t)| = 


-dy + CS [\^VMn^-^+4^Vn\^dT 
Jo 


Jz=o 4(l + |V^/i(i)P)® 

+ CsA^it) J^Ym + ||Vr;||^^_i + \\TZlj\\Y^dT + J g\V‘^{t)\'^dVu (4.2.30) 

and for Z°‘ = 9™, it holds that 

f Bdr <Ao[\ho\iim + cT\Vyho\lim + ||4^'"(0)f + |k(0)||^™_i] + Ao|Vy/i|^,„„i 
^0 

\VyZ^h{tr 




d‘^p{t)\Z‘^h{t)\ dy-a 


’z=0 


dy (4.2.31) 


.=0 4(l+|V,Mi)P)^ 

+ CS [ e^VMnm-i+e\\VV"^\\^dT + CsA^{t) [ + m\YmdT. 

Jo Jo 


Using (13.51) and integration by parts, one obtains that 


ft ft 

/ Jdr = 

'o Jo Js 

IQ“(0P 


Q“(-—- —V ■ V^Q“ + n^)dVrdT 

IP IP 

Js 2jp{t) * Js 2jp{0) ° 


< - 


IQ^(^)P 

2jp{t) 


■dVt+ f *^"^°^f dVo + Aoo(t) f ||g™(t)fdr+ [ / g“7^^dV.dr, (4.2.32) 

Js 27p(0) Jo jo js 


where 


Is 

l.„,o 1 . V 1 U 


7^^ = Z“r?(an—)arp+a^(—) • v» -c“(d) - —q(p) -c“(p) - [z“, —jsfp- [z“, —] • w 

7P 7p 7p 7p 7p 7p 


Furthermore, it is easy to obtain that 

f [ Q‘^Z‘^y{dt{—)dtp + dn—)-V^p)dYrdT<A^it) f WQ^itW + \h\j^^dT. 
Jo Js 'IP IP Jo 

It follows from (I3.7I) - (I3.8I) that 


(4.2.33) 


g“(-C“(d) - —C“(p) - — •C“(p))dV.dT 
lo JS IP IP 


<A(l,|HU,t)( / Iig^fdr 

rt 


i\\C-{d)\\ + \\C^{p)\\+C^{p)\\rdr 


<S [ ||Vp||^„_ldr + C5Aoo(t) f \\Vvf^rr.-l+\\Q^it)f + \h\nr.dT. 
Jo Jo 

Using (|2.2p . f|4.Q.5p and f|4.Q.6p . one gets that 

■ V'p)dV,dr 

< [' IIO”ll(ll[2”, -19711 + IIP”. —\ ■ V7II)* 

Jo IP IP 


(4.2.34) 


<Ao.(t)( / WQ-fdr] ( / (|b||^™ + ||Vp||^„_. + ||H|^™ + |/r|^.dr 
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[ \\'^p\\nm-idT+ CsA^{t) f \\v\\l^m + \\p\\l^m + WQ’^f + Ihll^mdr 
Jo Jo 

< S f ||Vp||^_,dr + CsA^{t) f ||F™(t)f + IIQ^Wf + \h\j^mdT. 

Jo Jo 


Combining (I4.2.32I) - (I4.2.35I) leads to that 

IQ“WI 


/o 


s ^ip{t) Js 27 p ( 0 ) 

t 


dVt-\-S [ WVpWl^m-idr 
Jo 


Therefore, using holder inequality, one gets that 


n‘^V°‘dVrdT 


0 Js 


< 


\\n-M\\\\V^\\dr 


and 


4 ' 2jp{t) 


dVt + l: [ + cr 

^ J z^O 


{i + \yyh{tW)i 


dy 


+ 11 2pe\SW°^\'^ + \e\div'^V^\^dVrdT 
Jo Js 


(4.2.35) 


+ CsA^{t) [ \\Vvf^^., + \\V^m^ + \\Q^it)f + \h\lir.dT. (4.2.36) 

Jo 

Now we estimate the terms involve TZ‘^. Using (13.611 - (13.1011 . it is easy to obtain that 

[ \\nljfdT<A^{t) [ e^\\V\f^^-^+Y^ + e\Z"^h\ldT. 

Jo Jo ^ 


<S [ e^\\V^v\\l^^-,dT + CsA^{t) [ + e\Z"^h\ldT, (4.2.37) 

Jo Jo ^ 

[ Wn'ljWY^dT <6 f e^\\V^vf^m-idT + CsA^{t) j Y^ + ||Vn||^™-i + e\Z^h\ldT. (4.2.38) 
Jo Jo Jo ^ 

Combining (14.2.411 . (14.2.3011 . (I4.2.36L (14. 2. 3711 and (14.2.3811 . we obtain, for Z“ ^ 

'2 1 /■ |V,Z“h(t)|2 


< Ao[|(/i,v^V,M(0)|?,™ + ||(U™,g™)(0)f + ||n(0)||?,™] +5 / ||Vp||^„_,dr 

Jo 

+ A 06 f eWVV'-^f+e^VMH^-^dT + A^{t) [ A{Y^) + \\Wv\\l^m-i + e\Z''^h\ldT. (4.2.39) 
Jo Jo ^ 

Therefore, summing over a, using the a priori assumptions (14.0.411 on the Taylor sign condition and 
Lemma [2.101 and taking S suitably small, one proves (14.2.111 for the case Z“ ^ 5™. 

Collecting (I4.2.4L (I4.2.3ip . (I4.2.36L (I4.2.37p . (I4.2.38L and using (14.2.111 for the case of Z“ ^ 
whose proof has already been closed in (I4.2.39L one yields immediately (14.2.Ill for the case Z“ = 9™. 
Therefore, the proof of Lemma 031 is completed. □ 


4.3 Estimates for and div'^v 

To deal with the compressibility, one needs to obtain some uniform estimates for the pressure and the 
divergence of the velocity field. For incompressible flows, the divergence of the velocity field is zero 
and the pressure plays a role of Lagrangian multiplier so that the pressure can be estimated by elliptic 
regularities, see m- While, the pressure p{q) = g^, shown in (11.3111 ^ for compressible flows, satisfies a 
transport equation due to (ll.31ll ^ so that it can be estimated by energy method. However, it seems that 
one can only expect uniform estimates of m — 2 order conormal derivatives of Vp, which is one order 
lower than the one in incompressible flows in [32] , since the divergence of the velocity held is not free but 
depends on the regularity of the free surfaces in compressible hows. 

Precisely, one can obtain the following estimates for V^p and div'^u. 
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Proposition 4.5 For any t G [0,r®] and m> 2, it holds that 


Jo 

< \\Vvo\\‘^m-2-\-Aoe\\{div'^^vo,\/'^°Po)\\l^m-i-\-{Si-\-SAo)s‘^ [ \\V‘^v\\‘^m~.idr 

Jo 

[ \\VV^fdT + CsA^it)O + \\^dtv‘^v\\l,0[ A(y™) + ||+ e|Z'"/i|idr, (4.3.1) 
0i Jo ’Jo ^ 


and 




[ llVpIl^^.idr < Aoe^ / WVMn^^^dr + A^{t) f A{Ym)dT, (4.3.2) 

Jo Jo Jo 


10 Jo 

where 5i > 0 and (5 > 0 are small constants which will be chosen later. 

It should be noted that the estimate of f* \\V^ Z"^~^d\v'^v\\‘^dT is the key to control || V^z;||^^_idT. 

This is not necessary for incompressible case, since f* || V^u||^„_iC?t does not appear due to (ll.62|) . 

The proof of this proposition is a consequence of the following lemmas. Note that 


[ l|Vu(t)||2^™_2 < \\Vvo\\'^^-2 + /o ||9tVu||^„_2dr < || Vuo||^„-2 + /J || V?;||^„_idT, 
\l|Vp(t)||^™-2 < ||Vpo||«™-2+/o ||atVp||^„_2dT < ||Vpo||^^_2+/o ||Vp||^„_idT. 


(4.3.3) 


(4.3.4) 


We will estimate elKV^p, div‘^n)||^m-i, /g IjV^Z™ ^div‘^u||^c?T and || Vp||^^_iC?T in the remaining 
part of this subsection. The estimate of WA/vW^m-idr will be left to later subsections. 

The first lemma gives the equations satisfied by Z°‘v and Z°‘W^p. 

Lemma 4.6 For |a| < m, the equations for Z°‘'V‘f‘p and Z°‘div^v read as 
IZ^div^v +j^{df + v V^)Z<^p = C^{p), 

\g{df + V ■ V^)Z“u + - Z°‘{2pediv^S‘^v + XeV^div^v) = C^(p), 

where the commutators are given by 

C?(p) := - — ([Z“, Vy]VyP + [Z“, v,]d,p + K[Z“, d,]p) - , —](9fp + n • V», 

IP IP 

Chip) ■■= -g{[Z^,Vy]VyV + [Z“, v,]d,v + P,[Z“, d,]v) - [Z“, g]{dfv + v ■ V^v). 

Moreover, these commutators vanish when |a| =0. 

Proof. The mass equation (11.311) ^ yields that 

div'^z; + —{dfp + v ■ V^p) = 0 . 

IP 

Applying Z“ to (14.3.51) and (ll.31|) g shows the lemma by trivial calculations. 

Lemma 4.7 For any t G [0,r®], m > 1, it holds that 

e||(c?iz;‘^?;, V‘^p)(t)||^m_i + f || 

Jo 

< Aoe\\{div'^°vo,^'^°Po)\\H^-i+h [ e^||V\||^„-i + (5 / ||Vp||^m-idr 

^0 Jo 

+ ^e [ ||VP™fdT + QAoo(t)£" [ ||V 2 z;||^„,_ 2 dr 
0l Jo Jo 

FC,,A^{t)ll + \\Vdiv’^v\\l^) [ A(y^(T))+£|Z’"/i|idT 

Jo ^ 

where 5 and Si will be chosen later. 


(4.3.5) 

□ 


(4.3.6) 
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Proof. Multiplying (14.3.4|) ^ by Z°‘Aw‘^v and integrating over S x [0,t] lead to that 

e [ [ g(df + V ■ • V'^Z“div‘^'i;dV^dT + e f f Z°‘V^p ■ Z°^Abf'^vdVrdT 

Jo Js Jo Js 


>0 Js 

= {2p + A)e^ [ [ Z°‘V'^diY^v ■ V‘^Z°^dw^vdVrdT 

Jo Js 

rt 


+ £ 


[ f Z“(V‘^ X (V^ X i;)) • V‘^Z“div‘^i;dV^dT 

Jo Js 

: [ [ Chip) ■ \/‘f’Z°‘div^vdVrdT = V J,. 

Jo Js 

Integration by parts shows that 

e [ [ gidf + V ■ V'^)Z°‘v ■ V^Z^div^udV^dr 

Jo Js 

: [ [ g{dt + V ■ V‘f’)dW‘^Z°‘vZ°‘dW‘^vdVrdT 

Jo Js 

- [ [ (V^p(5t +vyVy + V,d,)Z^v + 0(V‘^z;)*V^Z“i;)Z“div‘^i;dV^dr 

Jo Js 

rt 


= —e 


— e 


+ £ / / g{dt +Vy ■'^y + Vzdz)Z°‘v ■ NZ^^div'^tidVT-dr = Ki. 

'o Jz=0 1 


Note that 


Z[9r,Z“]/ = a,]Z“/ - Z(|^)a,Z“/ + [df, ZZ“]/, [Z3,5,] = - 


d 


It follows from (|2.1L (12.21) and (12.131) that 

Ki = -is J g|Z“div‘^n(t)pdVt + j 6 'o|- 2 “div'^'’u( 0 )pdVo 

— £ [ [ g{dt + Vy-Vy + Vzdz)[div‘^, Z°‘]vZ°‘div‘^'^vdVTdT 

Jo Js 

< -is [ glZ^^diY^viOl^dVt + is / ^.o|^“cliv‘^«^;(0)pdVo 

^ Js ^ Js 

+ Aoo(t)s [ ||(cit + r;j; • Vy + -^i^Z3)[div‘^,Z“]ti||||Z“div‘^r;||dT 

Jo ^ 

< -is£ ^?|Z“div‘^^;(^)| 2 dVt + is^ ^.o|^“cliv‘^%( 0 )pdVo 

+ Aoo(t) [ £||V?;||^„_1 + £|Z'"h||dr, 

Jo ^ 

where one has the fact I4|z=o = 0 and 

II—lloo.t < llP.lloo.t + ll^.y.iioo.t < A(l, ||(u,a,u)iu.t + 

Z Co 


Similarly, one gets that 


f Ikllw'" + + s|Z™ ^h||dr, 

Jo ^ 


(4.3.7) 


(4.3.8) 


(4.3.9) 


(4.3.10) 


K2 < Aoo{t) 


(4.3.11) 
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As to the boundary term K 3 , one can obtain from the trace estimates (12.51) that 
Ar 3 <Aoo(t)e [ \ZZ°v\L2(R2)\Z°‘div‘^v\L2(R2)dT 

Jo 


< Ac 


(||VZZ“z;||^||ZZ“Hl^ + \\ 2 Z^v\\) 

■ (||VZ“div‘^i;||^||Z“div%||5 + ||Z“div‘"i;||)dr 


< 6 ie^ 


\\V^v\\^^-^dT + Sie / WVV^rdT 

Jo Jo 

+ C’(5iAoo(t) f ll'cll^m + IIVu|||^m-i + e|Z™/ipdr. 

Jo 

where (14.3.101) and Corollarv l2.6l have been used. 

Substituting (14.3.91) . (14.3.111) and (14. 3. 121) into (I4.3.8L one has that 

e [ f g{df + V ■ V‘^)Z°‘v ■ V^Z“div‘^i;dVrdr 

Jo Js 

<--^e [ g\Z°div'^v{t)\^dVt +1:6 f go\Z°‘div‘^°v{0)\'^dVo + Sie^ f ||V^z;||^m_idr 
J Js J Js Jo 

+ Sie [ \\VV^fdT + Cs,A^{t) f ||u||^^ + ||Vz;||^„_i + |Z™/i|2dr. 

Jo Jo 

For the pressure term, it follows from f|4.3.5p that 


/o JS 
rt 


= e 


= e 


Z^V^p ■ V^Z^dW^vdVrdr 

: [ [ Z^V^p‘Z^V^div^vdVrdr^e [ [ Z^V^p[V^, Z^]div^vdVrdT 
Jo Js Jo Js 

■ [ [ Z^V^p- (N‘^,Z‘^]div'^v + Z°‘(V‘^(—)dfp + V^(—)-V‘^p)]dVrdT 

Jo Js ^ JP -fp J 

-el [ Z°‘V‘^p ■ Z^^i — idf + V ■ V‘^W‘^p)dVrdT =: Li + L 2 . 

Jo Js IP 

Using m and (12.131) . one obtains immediately that 

llWllw™-ic^T +Aoo(t)(l + ||Vdiv‘^v||^ t) / A{Ym)dT 

Jo Jo 

C5Aoo{t)£‘^ [ \\V^v\\l^m-2dr 
Jo 


Li < 


As to L 2 , note that 


+ V ■ V‘^)V» = —(df + V ■ V‘^)Z“V^p + —IZ^^, ddv^p 
7P 7p 7p 

+ y Jz^(F)z"5(V^p + Z^(^)Z"VyV> + Z^(—)Z"9,vy 

^ \ 7 P 'yj) •759 / 


/3+u—cy 

I/3|>1 


7P 


IP 


then, it follows from integration by parts and (12.101) that 

-6 [ / Z“V 7 -+ V‘^)Z“V 7 dV^dr 
Jo Js IP 

: [ ^|Z“V 7 prfVt + e [ ^iZ^V'^^oPdVo 
Js ^IP Js 274*0 


= —e 


(4.3.12) 


(4.3.13) 


(4.3.14) 


(4.3.15) 


(4.3.16) 
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+ e I I (5n^)+div‘^(^))|Z“V‘^ppdVtdr 


/o Js 

,| 2 , 


<-£ / _|Z“V‘^p|"dVt + e 

Js 27 P 


f ^iZ^V'^'^popdVo + AooW re||Vp||^™_i + 

^7Po Jo 




(4.3.17) 


Since 14 vanishes on the boundary, one has that 

rt 


— e 


[ f [Z“,5,]V>dV,dr< AooW [ e||Vp|| 

Jo Js 'IP Jo 


-1+e|Z™-^/i|ldT. (4.3.18) 


For the term like 


—£ 


c;3,,Z“V‘^pZ^(—)Z"5tV‘^p, 


1 


/o 7P 

where /3 and ly satisfy that /3 ^ 0 and /3 + ly = a, it follows from (EH) that 


— € 


: [ fc 0 ,,Z‘^V^pZ^{ — )Z’'ZoV^p<[ e\\V^p\\^m-idT + e f ||Z^-iZ(7)z"ZoV>fdr 

Jo Js IP Jo Jo IP 

<Aoo(t) [ e||Vp||^^„i + ||p||^„_i + £|Z™“^/i|idr. (4.3.19) 

Jo ^ 

In a similar way, one gets that 

: [ /c;3..^“V‘^pZ'5(^)Z"V^V^p< Aoo(t) f Ym{T)dT. 

Jo Js IP Jo 


— e 


To estimate the term like 


one can rewrite it as 


(4.3.20) 


1^^ 


-£ / / C;3,,Z“V‘"pZ^(—)Z"5,V‘"p, 


0 JS 


IP 


-e I ! c^_-Z“V>Z'^(d—4^)Z3Z"V>, 


0 Js 


jpz 


where |/3| + li/l < m — 1, 1/ < m — 2 and - are smooth bounded functions depend only on z. Indeed, 
this follows from the fact that Z3((l — z)/z) = c(l — z)/z for some smooth bounded function c. 

If /3 = 0, |£'| < m — 2, then (14.3.101) implies that 


— e 


' Z“V‘^p(7-^4a)Z 3Z"V‘^P< Aoo(t) [ e||Vp||^„_i+e|Z™-i/i|idr. 


/o JS IPZ 

When /3 7^ 0, it follows from (12.11) that 


(4.3.21) 


10 


— e 


0 JS 


Z“ V>Z^ (——dii )Z3Z^V^p 

7P2 

{ ||Z“V>f dr}4 ||Z^(4^^)Z3Z"Wf dr} ' 


< £ 


(4.3.22) 


where one has used (I3.22p . (13.231) and (14.3.lOF Using (12.11) and the Hardy inequality (Lemma 8.4 in [35]), 
one gets that 


r II < A^(t) r ||p||^_, + 

Jo iP^ Jo ^ 


< Aooit) f Wpfu^-. + \\VA\l^-. + ||VU.|||,™_idr. (4.3.23) 

Jo 
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It follows from (I3.16|) that 

/ l|VF,||^_ldT< Aoo(t) / \\v\\\^-,+\\Vvf^rr.-l+\Z^h\ldT. 

Jn Jo ^ 


Combining (I4.3.21I) - (I4.3.24I) leads to that 


—£ 


[ Aoo(t) f Y^{T)+e\Z"^h\ldT. 

Jo Js IPZ Jo " 


(4.3.24) 


(4.3.25) 


(4.3.26) 


Substituting (I4.3.15L (j4.3.17l) - (|4.3.2()l) and (I4.3.25|) into (14.3.1411 yields that 

e [ [ Z°‘V‘^pV‘^ Z°‘div’^vdVrdT 
Jo Js 

<-e f ^iZ’^V^pfdVt+e f ^\Z‘^W^°pofdVo + CsA^{t)e^ f \\V\f^m-2dT 
Js ^IP Js ^IPo Jo 

+ S [ ||V>||^™-idT +Aoo(t)(l + ||Vdiv%||^_J / Y^{T)dT. 

Jo Jo 

Now, we estimate the RHS of (|4.3.7p . It follows from (|2.2p (|3.22p and (|3.25p that 

Ji = {2fi^Xy [ [ iV^^Z^dW^vl^dVrdr ^ {2^^ X)e^ [ [ [Z^ ,V^]diY^v ■ Z^dW^vdVr dr 

Jo Js Jo Js 

■ y |V‘^Z“div‘^z;|2dV,dr-Aoo(t)(l +||Vdiv‘"7;||^_J \Z^-^h\\dT 

- Koo{t)e^ [ \\^‘^v\\l^m-2dr, 

Jo 


2a -f A n 
> —e 


(4.3.27) 


and 


iJsl < l{2p + X)e^ j\\V^Z^-^dW^vf dr + C j\\CUp)f dr 

<t(2/x + A)e^ [ ||VZ™“Miv‘^a||^dT +Aoo(t) / Ym{T)dT. 
o Jo Jo 

Using integration by parts, one can obtain from (12.1211 and the trace estimates (12.511 that 
IJ2I = Ate^ / [ Z“(V‘^ X (V"^ X a)) • V‘^Z“div‘^vdVrdr 

Jo Js 

<i(2/x + A)e^ [ ||V'^Z“div‘^'y|pdr + Ao [ e^||V^a||^„_2 + Aoo(t) [ Ym{T)dT 
0 Jo Jo Jo 

+ ps'^ [ [ (V*^ X Z“(V'^ X v)) ■ V^Z^dW'^vdVrdT 

Jo Js 

<I(2/r + A)e^ [ ||V‘^Z“div‘^?;||^dT + AqE^ [ ||V^?;||^„_2 + Aoo(t) [ Fm(r)dr 

o Jo Jo Jo 

+ pe^ [ [ (N X Z“(V^ X v)) ■ {Vy, d)Z°‘ddv^vdydT 

Jo jz^O 

<t(2/x + A)e^ [ ||V'^Z“div‘^'y|pdr + Ao [ e^|| V^a||^„_2 + Aoo(t) [ Ym(T)dT 

o JO ^0 JO 

+ [ |(N X Z“(V^ X v))\i ■ |z“div‘^a|idr 

Jo ^ ^ 

<]{2p + X)e^ [ ||V^Z™-^div%fdr + 5i£2 / ||V^vH^^-idr + / ||VU™fdr 

4 Jo Jo oi Jq 

+ Cs,Ao,{t) [ y^(T) + £lV2r;||^™_2+£|Z'"h|idr. 

Jo ^ 


(4.3.28) 


(4.3.29) 
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Combining (I4.3.7I1 . (I4.3.13I1 . (I4.3.26I) - (I4.3.29I) and summing over a(\a\ < m — 1) complete the proof of 
this lemma. □ 


In order to close the estimate in Lemma HTTl it remains to bound || We first estimate 

II VZ'"“^div‘^u|pc?T by using of the mass equation (14.3.41) . It should be mentioned that, in general, it is 
difficult to derive an uniform estimate for \\Z'^~‘^dzz'u\\‘^dT due to the possible appearance of boundary 
layers. However, one can expect the estimate of || VZ'"“^div'^w|pdT due to weaker boundary layer 
for divu. Moreover, this estimate is very helpful to bound e || V^n||^m„ 2 dT, which will be used later. 
Precisely, we can prove the following lemma. 

Lemma 4.8 For m > 2, it holds that 

[ ||VZ™-2d*^‘^^;||2dr< Aoo(t) f ||Vp(t)||^„_i+ y^(T)dT, (4.3.30) 

Jo Jo 


and 


■[ \\VMn^--dT + e^ [ \\V^V^vf^^..dT<A^{t) [ Ym{T)dT. (4.3.31) 

Jo Jo Jo 


Proof. For |Qi|<m — 2,m>2, it follows from (j4.3.4p that 

VZ“div‘^n = -V(F(at +vyVy + Vzdz)zy) 

- V{^i[Z^,Vy]VyP) - V([Z“, ^]idtP + Vy • VyP)) 

1 1 3 

-V{ — i[Z^,Vz]dzP + Vz[Z‘^,dz]p) + [Z^,—]iVzdzp)) =:y^J^■ 
IP IP ^ 

It follows immediately from (EH)-( 1121 ) and the fact 14 U=o = 0 that 


[ / |JlpdV.dT< Aoo(t) / ||(p,Vp)||2,_i, 

Jo Js Jo 


and 


lo JS 

To estimate J 3 , one first notes that 


V{[Z^,Vz]dzp) = cpA^zf^VzZ'^dzP + Z^VzVZ'^dzp) 


I/3|>1 

Due to EH and (14.3.241) . one can obtain that 


\\VZf^VzZ'^dzpfdT< Y. / II^^V14Z"5,pfdT 


\0\<m—2 


< Aoo(t) [ ||VI4||^m-2 -b ||9^p||^m-2dr 
Jo 

<Aoo(i) / \\v\\‘^,n-2 + \\Viv,p)\\l^^.2 + \h\l^mdT. 
Jo 


On the other hand, it follows from (12.11) . Hardy inequality, (13.251) and (I4.3.24|) that 


(4.3.32) 


(4.3.33) 


f /" |J2pdVrdr < Aoo(t) / ||(p,n)||^™_2-b ||(Vp, Vu)||^^_2dT. (4.3.34) 

Jo Js Jo 


(4.3.35) 


[ \\Z^VzVZ''dzpfdT< [ \\ZP{^-^Vz)ZZ‘'dzp\\^dT 

Jo Jo ^ 

< Aoo(t) [ II-I4||^m-2 -b ll^zPlI^m-ldr < Aoo(t) / 11 (Pz , V Pj, ) 11 ^m_2 “b 11 9^p| | 1 cIt 

Jo z Jo 

<Aoo(t) [ ||(n,Vn)||^™-2-b ||Vp||^™-i-b |/i|«mdr. 

^0 


(4.3.36) 
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Thus, 


< A„ 


< 


[ [ \V{ — [Z‘^,V,]d,p)fdVrdT<A^{t)[ [ |[Z“,y,]9,p|2 + |V([Z“,y,]a,p)|2dV.dr 
Jo Js IP Jo Js 

M [ \\V,fnm-2 + \\dM\nm-3 + \\^Z^V,Z’'dM\‘^ + \\Z^V,VZ''d,p\\^dT 

Jo 

AooW / \\{v,Vv)\\l^m-2+\\Vp\\l^m-l + \h\l^^dT. 

Jo 

Similarly, 

[ / |V([Z“, —])(14a^p)pdVT-dr < Aoo(i) / ||(u,p, Vu, Vp)||^™-2 + 

Jo Js 'IP Jo 

Finally, it follows from (12.2L p.23l) and (14.3.101) that 

[ / |V( —14[Z“,5,]p)|2dV.dr < Aoo(t) / ||Vp||^™_2dr. 

Jo Js IP Jo 

Combining (14.3.371) . (I4.3.38|) and (I4.3.39L one gets that 

f [ iJol'^dVrdT < Aoo{t) [ ||(u, Vu)||^™_2 + ||(p, Vp)||^„-1 + |/i|^™dr. 

Jo Js Jo 

Then, (I4.3.3n|) follows from (I4.3.32I) - (I4.3.34I) and (|4.3.40|) . 

By p.31|) o. (I4.3.30|) and Corollary 12.61 one obtains that 

[ \\d'^v\\l^m-2dT < f ||VZ'"“^div‘^u||dr +Aoo(t) [ Ym{T)dT 
Jo Jo Jo 

< Aoo(t) [ Yjn{T)dT. 

Jo 

Then, (14.3.311) follows immediately. Therefore, the lemma is proved. 

Now we can estimate || Vp||^„_idr in the following lemma. 

Lemma 4.9 For every m > 1, it holds that 

[ WWpf^m^idT < Aoe^ [ \\vMnrn~^dT + A^it) f A(r™(r)) + e|Z™- 
Jo Jo Jo ^ 

Proof. For any |Qf| < m — 1, it follows from (I4.3.4I) „. (12.IF (12.21) . and (12.121) that 

[ ||Z“V‘"pfdT< f Wgidt+Vy-Vy + V^dOZ^^vfdr 
Jo Jo 

+ £^ f ||div'^S"^u||^™-i + ||V'^div‘^u||^m-idr + f \\[Z°', g]{dtv + 

Jo Jo 

+ [ \\g{[Z‘^,Vy]WyV+[Z<^,V,]d,v + V4Z‘^,d,]v)fdT 
Jo 

< Aoo(t) f Wvfn^ + IIVu ||^„„2 + \h\l,m + Ib||^™dr + f ||V^V^u||?,„_.dT 
Jo Jo 

<Aoe^ [ ||V\||^™_idr + Aoo(t) [ Y^ir) + e\Z"^-^Vyh\ldT, 

Jo Jo ^ 


(4.3.37) 


(4.3.38) 


(4.3.39) 


(4.3.40) 


(4.3.41) 

□ 


(4.3.42) 


Vy ■ A/yV + Vzdzv)\\’^dT 


which combining (12.111) yields the proof of this lemma. 


(4.3.43) 

□ 


Since div'^u is not free, the following estimate is the key to estimate f* Ib^f H^m-idr: 
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Lemma 4.10 For m > 6, it holds that 


\\div‘^v\\l^,r.-i < A( —, Qm{t) + Q{t)), 
Co 

where Q„i{t) = \\{v,p){t)\\‘^m + ||(Vi;,Vp)(t )||^„_2 + 

Proof. Applying Z“ with |a| < m — 1 to (14.3.51) gives that 


Note that 


then it follows that 


z-dw^v = -z“(^ + ScZi-E) _ z«(h5f£). 

7P 7p 7p 


^ ifco P 




l/3|-[^] .. |/3|=m-l 

|/3|=0 ^ l/3| = l+[f] ^ 

|/3| + |j/|<m—1 l/3| + |i'|<m—1 


(4.3.44) 


(4.3.45) 


(4.3.46) 


(4.3.47) 


For the L°° norm, it holds, for |/3| < [^1, that 

E E l|V27l|i||^7l|2<|b||^[^,+. + ||Vp||^i™,+,. (4.3.48) 

|/3|<[f] |/3|<[f] 

Thus, for m > 6, that is [^\ + 2 < m and [y] + 1 < m — 2, one gets that 



m<m 

(4.3.49) 

which implies that 

ll^llE-i < A(-,Q„(t) + Q(t)), m>6. 

P Co 

(4.3.50) 

Similarly, one has that 

E \\Z^v\\l^<H — ,Qm{t) + Q{t)), m>6, 

rn 

(4.3.51) 


which yields that, for m > 6, 


l/3| = [f] 






l/3|=0 

|/3| + |!y|<m-l 


l/3| = l+[f] 

|/3| + |i/|<m-l 


< H-,Qm{t) + Q{t)). 
Co 


(4.3.52) 


Additional care is needed to estimate since it involves dzP- Rewrite this term as 

2<^(Y^) = Ylz<^dzP + dzpZ<^{^)+ E CkZ^{ — )Z‘'dzp. (4.3.53) 

P P P P 

0-\-v—a 

Then, it follows from (13.251) that 

W-Z-dzPr < ll^^^^llLlbllE < A(-, Q{t) + Q™(t)), (4.3.54) 

P zp ^ Co 
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and 


||5,pZ“(^)(i)f < A{-, Qit) + Qrnm 

P Co 


(4.3.55) 


\l3\—m — 2 


,Vz 


l/3| = l ^ 

^+1^=0; 

Y V \0\—m-3 

< ll^(f )llLl|5.p||?,™-2 + ||z5,p||i.||^||?,„_.+ ^ \\zP{-^)\\l\\z’'dM(^ 

P P |/3|=2 P 

0-\-v>—a. 

< Qm{t) + Q{t)) + \\^\\Y-^,^\\d.p\\Y-^ 

^0 P 

< A{-,Qm{t) + Q{t)). 

Co 


(4.3.56) 


Collecting the inequalities (I4.3.45p . (j4.3.50L (14. 3. 52^ - 114. 3. 56p shows (14.3.441) . Therefore, the proof of this 
lemma is completed. □ 


4.4 Normal derivative estimates: Part I 

In this subsection, we will focus on the estimates of normal derivative of v. To close the a priori esti¬ 
mates, one needs to bound ||9zu||^m_idT, e ||3^u||^™_2(iT and Jq in this and next 

subsection. Since 

dzV = (dzV ■ n)n + n(dzv), (4-4.1) 

it suffices to estimate dzV ■ n and n(dzv). First, one can bound the normal component of dzV in terms of 
div'^u as follows: 

Lemma 4.11 For any k G N, it holds that 

f llclzU ■ n||^fc < Ao f \\div'^vWh^k + dr + Aooit) f \Z’"h\\dT (4.4.2) 

Jo Jo Jo ^ 

Proof. It follows from the definition of df that 


div"^?; = dfvi + 92 U2 -I- dtvs = -;^dzV ■ N -|- 9iui -|- 92^2. 


which yields that 


dzV ■ n = 


dzP 


VT+W^ 

Hence, (14.4.21) follows from (12.11) . Lemmaand (14.0.51) . 


dzF 


(div"^?; — 9iUi — 92^2). 


(4.4.3) 

□ 


Next, we will estimate the tangential components of dzV. Note that 

V'^uN = —diifdiV — d2>pd2V + d'^v = — diipdiv — 82^820, 


8 zP 


and 


(V'^u)*N= (9fu-N,92%-N,9>-N)* = (9iz; • N, 92U • N, 0)‘+ — (9^u • N)N. 

OzP 

One can obtain immediately that 


n(9,u) = 


8 zP 


l + \VvP\- 


r(n(V‘^uN) n(9i¥>9iu -h 82 pd 2 v)) 


8 zip 


l + \VvP? 

8 zP 


l + \VvP\- 


r(2n(S'^uN) - n((V^u)*N) T\{8ip8iv + 82 P 82 V)) 
r(2n(S'^uN) - n{(9iu • N, 82 V ■ N, 0)*} -h Ii{ 8 ip 8 iv + 92V592u)). 


(4.4.4) 

(4.4.5) 


(4.4.6) 


As a consequence of 8 zV given in (14.4.31) and (I4.4.6L one have the following estimates for 8 zV: 
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Lemma 4.12 For any k it holds that 

r < Ao f \\{dtv^v,S^)f^, + \\V^+^\\'^dr + f \Z>^h\ldT (4.4.7) 

Jo Jo Jo ^ 

where := 11(5"^i;N). 

Furthermore, it follows from the expression of dzzV and (14.3.301) that 
Lemma 4.13 For any m> 1, it holds that 

e[ \\^Mr)\\n^-2dT < Aoe [ \\V^Z^-^Sn{T)fdT + A^it) [ r™(r)dr, (4.4.8) 

Jo Jo Jo 


and 


[ ||V2u(T)||^„_idT < Aoe^ / ||VF™)(T)f dr 

^0 Jo 

+ Aoo(t) [ Y^{T)+e\Z^-^VyhiT)\ldT. 

Jo ^ 


(4.4.9) 


The key is to estimate S^. Now, we follow the argument in |32| for the convention-diffusion equation 
solved by to derive the estimates of e || V‘^Z'"“^S'„(r)|pdT and J* || V'^Z™“^S'n(T)|pc?r. The 
only difference is that div'^’u will be involved but has been estimated in the last subsection. It should be 
noted that Sn satisfies the homogeneous Dirichlet boundary condition 


S'n = 0, on z = 0. (4.4.10) 

which follows from (11.331) and will be important for the following analysis. 

Now, we start with the estimate of S'n in 77"*“^ which bound e||V‘^Z'"“^Sn(T)|p(ir. 

Lemma 4.14 For any t £ [0,r'^] and m>2, it holds that 


\\S^ml^-.+e f \\V^Z-^-^S^{r)fdr 
Jo 

< Ao||Z“Sn(0)f-pdAoe^ [ \\{V‘^Z^-^div^v,V‘^Z^-^Sn,V^V"^){T)\\'^dT 

Jo 

+ CsA{-,Q{t)+e\\A/^div'^v\\'L.,^) [ Y^{t) + e\Z^h\\dT. (4.4.11) 

Co Jo 

Proof. We first derive the equations satisfied by Sn. Applying to (ll.31|) 2 gives that 
gd'fV^v + q{v ■ ■ V^v - g,e/A^(y^v) = (/r -f A)e(V^)2div‘^u - {V^fp - g ® dfv - V^vV^[gv), 

where f denotes the Hessian matrix of /. Thus, 

pans^u) -b g{v ■ V^)(S^u) - pe/A^^S^v) = (m + A)£(V'^)2div^r; - (V^)^^ 

- ® ® V‘^p) - ^(V^vV^igv) + {V^vV^{gv)y )- 


Therefore, Sn satisfies 
where F is given by 
with 


gdfSr, + giv • V‘")Sn - peA^S^ = F, 
F = Fi + F 2 + F^, 


(4.4.12) 

(4.4.13) 


Fi = g{dfU -b V ■ V^n)(S‘^'i;N) -b pn(S‘^u(afN -b v ■ V^N)) (4.4.14) 

- in{(v^p(g)aft>-bafu(8) v‘^p)n} - in{(v‘^uV'^(H + (v'^uV'^(ez;))*)N|, 

F 2 = -n{(Vy, 0)*(V> • N) -b (V^N)*V>}, (4.4.15) 

Fs = (^-b A)£n{(Vy,0)*(V'^div‘^?;-N) - (V'^A7)*V^div‘^u} 

- At£(A^n)(S^nN) - 2/ieafnaf (S^uN) - ^enjs'^uA'^N -b 2afS'^uSfN}. (4.4.16) 
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where the summation convention has been used in the above expressions. It follows from Lemma 12.1 
Lemma [231 and (14.3.3111 that 


[ ||Fi||«™-2dT < Aoo(t) / ||(V?;,Vp)(r)||^™_2 + ||(p,u)||^„_i + |Z”" ^h\ldT, 
Jo Jo ^ 

f < K^{t) f ||Vp||?,._, + |Z™-i/r|idT, 

Jo Jo ^ 

f mw^.^dr < Ao f \\v\\l^.. + eWVMlirn-. + e^||V\||^„_.dr 
Jo Jo 

+ A(-,Q(t) + e||Vdiv %||2 f\Z-^-^h\\ +e\Z^h\\dT. 

Co do 


(4.4.17) 

(4.4.18) 

(4.4.19) 


The lemma is proved by induction. For the case m = 2, multiplying (14.4.1211 by S-^ and then using 
Lemma [231 and the boundary condition (I4.4.6p . one gets that 

f ^g{t)\Snit)\^dVt + ^ie f f \V^Sn\^dVrdT= f ^go\Snm^dVo+ f f F ■ S^dVrdr (4.4.20) 

It follows from (I4.4.17|I - (I4.4.19I1 with m = 2 that 

f F-SndVrdT< ( /'‘||(Fi,F 2 ,F 3 )fdT)'( /Vnfdr)' 

*-'0 'J 0 

< f + e^llVdiv^ull^, + e^ll 

Jo 

+ A( —,Q(t) + ||V^div‘^u||^_t) f |k||«2 + ||(p,Vu)||^i + |Z/i|i+e|z2/i|i 

Co Jo 2 2 

< A(-,Q(t) + ||V‘^div%||L_,) / A(y^(T))+e|z2/i|idT, (4.4.21) 

Co Jo ^ 

where (14.3.3111 has been used in the last inequality. Substituting (14.4.2111 into (I4.4.20L one gets from 
(14.3.3011 and (14. 4. 911 that 

£ igiOlSniOl'^dVt J \^‘^Sn\‘^dVrdT 

< f lgo\Snm^dVo+Ao£^ [ ||V‘"Zdiv^«,V‘"Z5„,VFidT 

^ S 0 

+ A(l,Q(t) + ||V^div‘"u||^_,) / A(y™(r))+e|Z2/i|idr. (4.4.22) 

Assume that (14. 4. lip has been proved for /c < m — 3, we shall prove it for /c = m — 2. First, applying Z“, 
with jal = m — 2, to (14.4.1211 yields that 


gdtZ°‘Sn + gv ■ V‘^Z°‘Sn - /ieA‘^Z“S'n = Z“F + Cs, 
with Cs = Cg + C| given by 


(4.4.23) 


= -[Z^,g]dtSn - [Z‘^,gvy]VySn - [Z“,eF,]5,5„ - (4-4.24) 


2=1 


and 

C| 

Multiplying (14.4.2311 by Z’^Sn with |q;| = m — 2 and integrating give that 

ri 

L2 


(4.4.25) 


git)\Z‘^Snit)\^dVt+fi£ [ [ \V^Z‘^Sn\'^dVrdT 

Jo Js 

j i^^o|^“^n(0)pdVo+ [ [ |Z“F.Z“5„|dV.dT+ f j Cs-Z^S^dVrdT. 
J S J 0 'J S J 0 'J s 


(4.4.26) 
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As a consequence of (I4.4.17ll - (|4.4.19ll . (I4.3.42I1 . (14.4.811 and (14.4.911 . we obtain that 


0 Js 


\Z^F-Z<^Sn\dVrdT<{ I ||(Fi,F2,F3)||^^_2dr 


)' 


\\SXu^-.dT 


< 5 f II+ e^\\y\\\nr.-^dT 

^0 

+ A(-,Q(t) + e||Vdiv^u||2 [ Y^{t) + e\Z^h\\dT. 

<^ [ \\V‘^ Z^-^Snf dT + S Aoe"^ f ||(V‘^Z™-Miv‘"w,V‘^Z”^-^5n,V‘^y™)(r)fdr 

4 Jo Jo 

+ QA(-,Q(t)+e||Vdiv%||2 f F„(r)+e|Z™h(t)|ldr. (4.4.27) 

Co Jo " 

It remains to estimate the terms involving Cg, C|. Similar to [32) . integrating by parts and using the 
Hardy inequality, one gets that 


[ [ |C^Z“5„|dVtdr < [ K{Yr^{T))dT. 
Jo Js Jo 


(4.4.28) 


(4.4.29) 


(4.4.30) 


For the term involving C|, we notice that A’^ can be rewritten as 

with the matrix E defined by 

This yields immediately that 
with 

Cf = Me[Z“, ^]V • (EWSn), Cf = /re^[Z“, V-](F;V5„), Cf = • [Z“, AV]^„. 

Oz^ 

By similar arguments as |32j . one obtains that 


A-/ = ^ 

—V • 

zP 

(EVf), 


f dzip 

0 

-dip 

E = 

0 

dzP 

-d 2 P 



-d 2 p 

I+IV^V’I 

dzV 


^2 _ /^ 2,1 I ^ 2,2 I ^ 2,3 


ClZ'^S^dVrdT\ < ^ / / |VZ’"-^,S„|"dV, + Ao£ 


0 JS 


VZ'"-^5'nrdV^dr 


+ Aoo(t) [ Yjn{T)dT. 

Jo 


(4.4.31) 


Plugging (j4.4.27L (j4.4.28ll and (14.4.3111 into (I4.4.26p . we obtain that 

j ]-Q{t)\Z^S^{t)\^dVt+^Jie j f |V‘"Z“5„|2dV.dT 
«/ ^ 0 ^ s 

< jT ipo|^“^n(0)pdVo + Ao£^ J \VZ^-^Sn\^dVrdT + A^{t) A(y^(T)) + £|Z'"/l||dT. 

Thus, the lemma is proved by (12.1011 and the induction assumption to control Ao£ || V‘^Z'"“^5'n|pd'r. 

□ 


The following lemma give the estimate of f* HV^Z™ ^S'nlPdT- 
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Lemma 4.15 For any t € [0,T®] and m> 1, it holds that 

eWS^mur.-, + f II dr 

< ^A(-,r™(0)) + Aodie2 [ [\V^Z"^-^dtv'^v\'^(WrdT 
Ol Co Jo J 

+ Ci.AooW [ A{Ymir)) + e\Z^h\l + ||Vz;||^„_idT, 

Jo ^ 

where di > 0 and d > 0 will be chosen later. 

Proof. Multiplying (14.4.231) by £Z°‘Sn with |q;| = m — 1 and integrating give that 
f yit)\Z^Snit)\^dVt+yi£^ f f |V‘"Z“5„|2dV.dr 

= [ ^go\Z^Sn{0)fdVo+£ f f Z^F-Z^SndVrdr +£ [ f Cs-Z^SMdr 
Js ^ Jo Js Jo Js 

Similar to (14.4.171) . it is easy to obtain that 

f ||Fi||?,„_.dr < A^{t) f ||(V^, Vp)||^„_. + ||(p,^)||^™ + |Z-h|idr, 

Jo Jo 


which implies that 


Z“Fi • Z°^S^dVrdT 


0 JS 


< Aoo(t) / Ymij) + £\Z'^h\\dT. 
Jo " 


Using integration by parts, one can get that 


Z°‘F2 ■ Z°‘S^dVrdT 


0 JS 


<£ [ ||VU™fdr + Aoo(t) [ Yrnir) + £\Z^h\ldT, 
Jo Jo ^ 


(4.4.32) 


(4.4.33) 


(4.4.34) 


(4.4.35) 


/o Js 

.2 /-t 


< 


and 


Z“F3 • Z^^SndVrdT 

^ lo I + / |V^Z™-Miv^?;pdV^dT 

^ f WVV^fdT + CsAooit) f\(Y^iT))+e\Z^h\ldT 

Ol Jo Jo ^ 

[ f ClZ‘^SndVrdT +e [ [ ClZ°‘SndVrdT 

Jo Js Jo Js 

J J\V^Z’^Sn\‘^dVrdT + Aa,it)J A(y^(T))+e|Z™h||dr. 


(4.4.36) 


/o JS 

Te 


(4.4.37) 


Substituting (j4.4.34l) - (l4.4.37l) into (j4.4.33p . we prove (I4.4.33p . Therefore, the proof of this lemma is 
completed. □ 

Then, suitable choice of S and di yield the following proposition 
Proposition 4.16 For any t G [0,r'^], to > 5, it holds that 

IKU"*, QmtW + \{h, v^V,h)(t)|2,„, ^ V^;)(t)||2,._. + eUdtv^v, V^p, 5„)(t)||?i™-i 

+ e|Z™h|i+ / ||Vp||^™_idr + e / ||Vc||^™dr+/ e||V^ull^™-^ + V^H^^-idr 

^ Jo Jo Jo 

< A(l,r™(0)) + Aoo(t)(l + ||V‘^d*t;%||L,t) [ A(r™(T))+e|Z™h|i + ||Vu||^„_idT. (4.4.38) 

Co Jo ^ 
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Proof. Combining Propositions 14.4114.51 Lemmas 14.131 in4l and l4.151 one can obtain that 


+ r£||VP™(T)|p+e||V\(r)||^_.+£2||V\(r)||^._.dr 
^0 

< j-A(l,y™(0))+Ao(Ji+Ao|-)£2 [ \\vMr)\\nm^idT + {S,+Ao^)e [ \\VV^\\dT 
0i Co di Jo Jo 

+ CsA^{t){l + ||V‘^div‘^n||ic» ) [ A{Ym{T)) + e\Z'^h\l + ||(4.4.39) 

Jo ^ 

Choosing (5i > 0 to be small such that and ^ ^ yields that 

||(P™, Q™)(t)f + \{h, y^Vyh)(t)\lr. + II Vn(t)||^„_. + e||(div‘^^;, V>, 5„)(t)||?,™-i 
+ f ell V‘^P-(r)f + ell V\(r)||^„_. + e^H V\(r)||^„_.dr (4.4.40) 

< A(-,r^(0)) + Aoo(<)(l + IIV‘"div‘"i;||2 ) [ A{Yrr,{T)) + e|Z™h|i + ||V?;||^^^iJr. 

This, together with (14.3.21) and Lemma [2^ yields immediately (I4.4.38L Thus, the proof of this proposi¬ 
tion is completed. □ 


4.5 Normal derivative estimates: Part II 

Now we start to estimate || Vt!||^m-i Jr. It follows from Lemmas 14.101 and 14.121 that one needs only to 
bound ll^nll^m-i. However, this seems difficult since HV^pH^m and ||V‘^div‘^n||?^m, which appear in F 2 
and F 3 in the equation (14.4.121) for ^n, cannot be estimated uniformly in e:(c.f. Lemma [4.1 5 1 ) . As in [32], 
we use the vorticity instead of Sn to perform this estimate. Let us set 

Wn :=a; X N = n(w X N) (4.5.1) 

where uj is the vorticity dehned by w = x v. Since (14.4.51) yields that 

2Sn = LOn + 2n{{div ■ N, d 2 V ■ N, 0)‘}, (4.5.2) 

one can obtain from (|4.4.10l) that 

(Wn)'' = -2n{(aiu • N, d2V ■ N, 0)‘}. (4.5.3) 


which implies that 

(Z“a; X N)'’ = -(w x Z“N)'’ - ([Z“, wx, N])^ - 2Z“(n{(ai?; • N, ^ 2 ^^ ■ N, 0)*}). (4.5.4) 


Thus, one can obtain the sharp estimate that on the boundary, for a < m — 1, 


y/e f |(Z“a; X N)'’pJr < Ao\/e [ |a;^|^m-2 -|- |c''|^m Jr-I-Aoo(t) [ v^|VyJr 
Jo Jo Jo 

<Aoe [ ||V\||„„.-2 + ||V?;||?i^+Aoo(t) [ r™(r)Jr. (4.5.5) 

Jo Jo 


Hence, one can only expect to bound HZ'" ^oj x N||^Jr by the similar argument in [32] ■ Here and 
henceforth, ||Z'"“^a; x N|p = X]|a|<m-i 11-2^“^ ^ N|p. It should be emphasized that, one may not 

expect to bound directly Z’^~^uj like m, since it involves y/e /J |a;^|^„_i Jr which seems impossible to 
be bounded due to the invalid of (11.621) in the compressible flows. However, as shown in (I4.5.2L Wn 
involves only the tangential derivatives of v so that one can follow the argument in |32| for Wn instead of 
w. Actually, we will bound x N to avoid too much regularity of h involving. 


The main result of this subsection is the following proposition; 
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Proposition 4.17 It holds that, for any m > 6 that 


where 


||Vr||^™-i)' <Ao||Z“wof+£ J + + ||V^Z™- 25 „fdr 

pt pt 

+ A{M{t)) / A{Ym{T))dT + A,^{t)ti / ||Vp||^™_idr. (4.5.6) 

Jo Jo 


M{t):= sup (A{-,QiT) + Qm{T)+e\Z"^h\\) + e [ ||VP®f + ||V\||^4ds'). (4.5.7) 

0<T<t ^ Co ^ Jo ^ 


The proof of this proposition is a consequence of (14.5.101) and Proposition 14.181 below. 
Since 


and 


Z“a;„ = Z“(w X N) = Z“w x N - [Z“, Nxjw, \a\ < m - 1, 


II - [Z“,Nx]o;f < ||u;||L||V,7y||?,„_, + ||c.||?,4,.o||V,r?||?,_. + ^ HZ'^N x Z^c.|p 

l/3| = l 


(4.5.8) 


(4.5.9) 


< A(-,Q^(t) + Q(t)), 

Co 

we can obtain from (|4.4.1I1 (I4.4.8|) . (14.4.61) (14.8.441) and (14.5.211 that 

[ I|i9zv||^™-idr < Ao [ ||Z'"“^w X N||^ + ||div‘^r;||^™-idr + Aoo(t) [ A{Ym{T))dT 

Jo Jo Jo 

<Ao [ WZ'-^-^ujxNfdr + A^it) [ A{Yrn{T))dT. (4.5.10) 

Jo Jo 

So, it remains to control Z“a; x N, for |a| < m — 1. 

Now we first derive the equation solved by Z°‘uj x N as follows. Applying V'^x to (11.311) ^ yields that 


gOfuj + gv ■ V'^uj — g,eAYuj = Fi, 


where F-\ 


Fi = —\J^g X dfv — V'^g x ((v ■ V‘^)v) + geo ■ V^v — gdiv'^vco. 
Applying Z“, with |a| = m — 1, to (14.5.111) gives that 

gdtZ'^uj + gvy ■ VyZ°‘uj + gV,d^Z°‘u} - g,eA'^Z°‘u} = Z°‘Fi + F2, 

where 


(4.5.11) 

(4.5.12) 

(4.5.13) 


F2 = [Z“,^^]atW+[Z“,^)^;y]Vya; + [Z“,ey,]5,a; + 0l4[Z“,a,]w + /re[Z“,A‘^]a;. (4.5.14) 

Multiply (14.5.13p to get that 

gdt{Z°^uj X N) + gVy ■ Wy{Z°‘uj x N) + gVA{Z°^oj x N) - g,eA'f’{Z°‘uj x N) = F, (4.5.15) 

where the source term F is given by 

F = Z“Fi X N + P 2 X N + P 3 , (4.5.16) 

with F3 given by 

F3 = gZ'^uj X 9tN + gVy ■ Z“w x Vj,N + gV^Z'^uj x d^N + neZ°‘uj x A'^N + 2/reV‘^Z“w x V‘^N. (4.5.17) 

Next, we follow the argument in [32], where Z™“^w is bounded, to perform the estimate of Z™“^w x N 
as follows. 
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Proposition 4.18 For any t G and \a\ = m — 1, there exists A{M{t)) such that 

llZ^io X N||i4(o,t;i.(5)) < Ao||Z“c.of + Aoe f || + W^Mn^-^dr 

Jo 

+ A(M{t)) f A{Ym{T))dT + Aao{t)t^ [ ||V(p,'!;)( 

Jo Jo 


^ -idr. 




(4.5.18) 


where M{t) is defined in (14.5.71) . 

The proof of this Proposition will be a consequence of the following lemmas. As in m. we split Z°‘u! X N 
into two parts of the form 

Z“c.xN = a;))+<^, (4.5.19) 

where solves the non-homogeneous problem: 


r gdt^Zh + g^v ■ + g^zOz^Zh - = F, 

I X No, (4.5.20) 

[<h\z=O = 0, 

while LoZ solves the homogeneous one: 


r gdtOjZ + gvy ■ VyUjZ + gVzdzOjZ - = 0, 

s w“|t=o = 0, 

X N)^ 


(4.5.21) 


The solution to (14.5.2011 can be estimated by standard energy estimates. 

Lemma 4.19 For any |a| < m — 1, it holds that 

+ /‘||V‘^<,fdr<Ao||Z“^of+ Ao£ f\\VMn^-zdr 

Jo Jo 

+ Aoo(t) / e\Z'^h\l + \\V{p,v)f^^-^+A{Y^{T))dT. (4.5.22) 

Jo ^ 

Proof. Noting the homogeneous Dirichlet boundary condition for one deduces from (14.5.201) by a 
standard energy estimate that 

i [ g\LoZh\^dVr + pe f [ |V^a;“JdV,dr < Aq [ ^,o|Z“a;opdV, + | f [ F ■ ojZ^dVrdr\, (4.5.23) 
^ Js Jo Js Js Jo Js 

where F is given by (|4.5.16p . Then (|2.H) implies that 

I f f Z<^F, X N . eoZndVrdrl + f j IF 3 • <H?dVrdT (4.5.24) 

Jo Js Jo Js 

-1^ j \^nhiT)\'^dVr+A^it) j ||(p,?;)||^™ + ||(Vp, Vu)||^™-1 + |Z’"-^/i|idT. 
tDo<T<tJ Jo ^ 

Using similar arguments as m, we can obtain that 

f [ \F2xN-cuZn\dVrdT 
Jo Js 

- t //+ j K^(r)|2dV. 

+ Aoo(t) [ e\Z"^h\l + IIV(p, u)||^„_i + AiYmiT))dT. (4.5.25) 

Jo ^ 

Substituting (14.5.2411 and (I4.5.25P into (14.5.231) proves Proposition 14.191 Thus, the proof is completed. □ 


To treat ujh, we will modify the approach of microlocal analysis used in m to obtain that 
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Lemma 4.20 For any \a\ = m — 1, it holds that 



\\VV^f + ||V\||^„- 2 dT + A{M{t)) 



A(Yjn{T))dT. 


(4.5.26) 


where M{t) is defined in (14.5.71) . 

In order to prove the above lemma, we follow the argument in [32) to use Lagrangian coordinates to 
eliminate the convection term so that it is convenient to perform the microlocal symmetrizer method. 
Let us define a parametrization by 


dtX{t,y,z) = u{t,X{t,y,z)) 


^oX), X{0,y,z) = <^>{0,y,z), 


(4.5.27) 


where <I>(t, •) ^ stands for the inverse of the map $(t, •) defined by (11.2411 . Define J{t, y, z) = \detXX{t, y, z)\ 
to be the Jacobian of the change of variable. The following estimates for X are proved in Lemma 10.5 in 

m- 


Lemma 4.21 It holds that for t G [0,T], 



||VX(t)||L=o + ||5tVX(t)||L=o < 

||VX(t)||i,oo + ||5tVX(t)||i,oo < A{M{t))e^XM{t)) ^ 

v^llV^Xlli.oo + V^\\dtX^X\\L<^ < A(M(t))(l + t 2 )e‘A(M(i)). 

(4.5.28) 

(4.5.29) 

(4.5.30) 

Now set 

0 X) 

(4.5.31) 

where 7 > 0 is 

a large parameter to be chosen. Then solves in S the equation 



ao(5tD“ + 7D“) — pLedi{aijdjTl°‘) = 0, 

(4.5.32) 

where 

oo = £>(t, oX)|J|5, (ay) = |J|5p-\ P,j=dtX-djX. 

(4.5.33) 


Thanks to Lemma [4.211 the equation (14.5.321) is a parabolic equation. On the boundary, it holds that 

■■= e-^‘(Z“w X N)(t, ($-1 o X){f y, 0). (4.5.34) 


The following theorem holds: 

Theorem 4.22 (Masmoudi-Rousset [32]) There exists 70 depending only on M > 0 such that for 
1 the solution of (14.5.3211 with the boundary condition ()4.5.34ll satisfies the estimate 

< A{M{T))V-ej^ (4.5.35) 

where the norm Hi {0,T; L^) is defined by 




= inf{||P/|| 


HZ{R,L^{S)) 


yPf = f on [0,r]x5}, 


with the norm on the whole space by Fourier transform in time. 

Proof of Theorem 14.221 The proof of this Theorem is almost the same as the proof of Theorem 
10.6 in [32] where the symmetrized method and paradifferential calculus are used. This is because the 
coefficients of (14.5.321) given in (14.5.3311 is the same type (determined by the estimate in Lemma 14.211 
II (i?! II 1,00 < M and initial data which ensure $0) of symbols as the one in [ 32 ] which is essential in 
the progress of taking use of paradifferential calculus. We refer to Theorem 10.6 in [32] for more details. 
□ 
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Proof of Lemma 14.201 As a consequence of Theorem 14.221 and the Sobolev embedding inequality, one 
has that 

rT 




2 

L4(0,T;L2(5)) 




i/4 (o,T;L2(5)) 


Jo 


Consequently, it follows from changing of variable and (14.5.51) that 


ll<lli^(0.t;L^(5)) < m{t))V~e f X 

<e f ||Vu||^^ + ||V2n|||,™-2dr + A(M(t)) [ Y^{T)dT. 

Jo Jo 

which yields Proposition 14.201 Therefore, the proof is completed. 


(4.5.36) 

□ 


4.6 L°°-estimates 


In order to close the estimates, one needs to bound the L°° norms of v, h and p contained in Aoo(t). We 
start with the following estimates through standard Sobolev embedding theorem and anisotropic Sobolev 


embedding theorem. 

Lemma 4.23 For any t £ [0,r'^], the following estimates hold: 

e N, (4.6.1) 

Up,vm\\n^,^<M-,Q5{t)), (4.6.2) 

Co 

\\Vp{t)rn.,^ < WA^pmh + M-, Q5(i))llVp(t)||?,4, (4.6.3) 

Co 

< ||AM0II«4 + A(l,Q5(t))(l + WVpitWn.), (4.6.4) 

Co 

WVdzv^vml < II+ A(l, Q,{t) + IIV^(t)||L)(l + IIVp||?,4), (4.6.5) 

Co 

\\Vdzv^v{t)f^^,^ < S\\A^p{t)f^. + CsAi-, Qeit) + ||Vn(t)||^i,oo)(l + ||Vp||?,5). (4.6.6) 

Co 

Proof. (14.6.11) follows from the two dimensional Sobolev embedding. The anisotropic Sobolev embedding 
( 1 ^ yields that 

\\{p,v){t)\\l.,^ < ||(Vp,Vn)(t)||„ 3 ||(p,u)(t )||„4 < A(l,Q 5 (i)). (4.6.7) 

Co 

To prove (14. 6. 3D . one notes that 


A'^ = 


l + |Vy¥>P 
\dM^ 


dzz + 


E (a- - 


2 di‘p 

dzP 


didz 


2 a(|^)& + + 4- 

OzP OzP OzP OzP 



which follows from the definition of df immediately. Thus, (12.4|) implies that 

l|Vp(t)||L < ||Vp(t)||3||9.Vp(t) + ||Vp(t)||i|| < ||Vp(t)||3(||a.p(t)||3 + WdzzPim) 

< A(l, |Mt)|2.oo)||Vp(t)||3(||Vp(t)||3 + Ib(t)||2 + II A>(t)||) (4.6.8) 

<||A>(t)f+A(l,Q5(t)). 

Co 

Similarly, 

\\2yp{t)\\l < ||ZVp||3||5.ZVp|| + llVpIll 

< A(l, |h(t)|3,oo)||Vp(t)||„4 (llVp(t)||„4 + Ib(t)||„3 + \\A^pit)\\nz) 

<||A>(t)||?,.+A(l, 25(0)11 Vp||?,4. 

Co 


(4.6.9) 
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Therefore, (14.6.311 holds true. Next, it follows from (14.3.511 . (14.6.211 and (14.6.311 that 


l|div%(t)||^.,=. < II 


< nil' 






< 


A(^, Q5W){i + ||Vp(t)||„4 [||Vp(t)||«4 + Ib(t)||„3 + ||A‘^p(t)||„4] } 


<||A‘^p(t)||^.+A(-,Q5(i))(l + ||Vp|||,4) 

Co 


Furthermore, it follows from the relation 


Vdiv'^u = - —(9tVp + Vy ■ \7yVp + 1492Vp) - V( —)9tp - V( —) • Vyp - V(—)dzp, 

ryp ryp ryp r^p 


that 


l|Vdiv‘^u(t)||^ 

< Ao{||Vp(t)||?, 4 ,. + lk,(t)||Ll|Vp(t)||?, 4 ,o. + ||(h^. 92 Vp)(t)||L + l|Vp(i)llLlbWII«u=» 
+ (llvu,(t)iiL + iiKvp)(t)iiL)ibb)ii?,i.=o + (iivy2b)iiL + ll(KVp)(t)iiL)ii92Pb)ii^ 

< A(l, Q5b))(l + II Vu||L)(l + II Vpllliu,.) 

Co 

< II A'^pII^, + A(l, Qsb) + II Vu(t)||^)(l + II Vp||?,4), 

Co 

since 14 vanishes on the boundary. Similarly, 

II Vdiv‘^u(t)||^4,„ < <511 AVb)ll«^ + Qeb) + II Vu(t)||?,4,o„)(l + || Vpll?,^). 

Co 

Therefore, the proof of this lemma is completed. 

We now estimate ||A‘^p||^i. Applying div"^ to the momentum equations yields that 
— {2p + A)eA'^div'^?; + A^p = div‘^{gdfv + gv ■ V‘^v) 

This, together with (14.3.511 and (13.1511 shows that 

edtA'^p + A'^p = div‘^(g9fu + gv ■ V^v) - e(vyVy + V^dz)A‘^p 

Zfl H" A Zf! A 

- e-fp(d^pA^{—) + V^pA^l — ) + 2V^9fpV‘^(—) + 2V‘^( —. 

V 7p 7p 7p 7p / 

Lemma 4.24 For m > 6; it holds that 

sup (||A>(r)||7+£||A>(r)||?,.)+ f WA^p^n^dr 

0<r<t ^ ''Jo 

< CAo(||(A»(0)||^i +e||(AV)(0)||«2) + Aoo{t) K(Ym{T))dT. 


0 

(4.6.10) 


(4.6.11) 

(4.6.12) 
□ 

(4.6.13) 

(4.6.14) 


(4.6.15) 
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Proof. Applying Z“(|q;| = 2) to (14.6.14L multiplying the resulting equation by and then 

integrating over S x one can get that 




[ [ Z^{-fpA^p)Z^A^pdVrdT 

Jo Js 


-e [ f Z‘^hpdtpA‘^(—))Z‘^A^pdVrdT - e [ f Z°‘(jpV^pA'^(—))Z°^A^pdVrdT 
Jo Js IP Jo Js IP 

-2e [ [ Z“(7p(V^9fpV^( —) + V‘^(—)V‘^V^p))Z“A>dV^dr 
Jo Js IP IP 

t 5 

■ [ / Z“{7pdiv‘^(^>(a,^+u-V‘^)u)}Z“A>dV.dr =: J||(Z“A‘^p)(0)f+ y J, 
Jo Js ^ 


(4.6.16) 


2p A 

The second term on the LHS above can be estimated directly as 
1 


2 . +a / Pp\Z~A-pfm4r 

+ —-— [ [ jZpZ{A‘^p)Z’^{A‘^p)dVrdT +—Z 

Jo Js 2 /r + A 


> 


2/i -(- A 
3 


■ f f 7 Z>A>Z“(A^p)dV^dr 
Jo Js 


> 


4(2/r + A) Jo Js 
3 


4(2/r + A) 


f [ 7P|2“A>pdV.dT-Ao /‘||Zp|lLl|A‘^p||?,A + ||Z>||L||A>fdr 
Jo Js Jo 

[ / 7p|Z“A‘^ppdVrdT - / A( —, Qm(t) + ||A‘^p||^i)dT, (4.6.17) 

Jo Js Jo Co 

where m > 5. Next, the terms on the RHS of (14.6.16|) can be estimated separately. First, it follows from 
integration by parts, (12.211 and (I3.25|) that 

I Ji| < £| / f {vyVyZ'^A^p + 14a,Z“A‘£p)Z“A>dV.dr| 

Jo Js 

+ £| [ [ {[Z^,VyVy]A‘^p+[Z^,V,]d,A^p + V,[Z^,d,]A^p)Z^A^pdVrdT\ 

Jo Js 

< A^{t) feWA'^pWl.dr + e f \\Zlz-V.\\U\A^p\\\, + \\ZlzVA\U\J^^p\\l.^dr 

Jo Jo ^ ^ 

< £^ / llV^ulll^adr + Aoo(t) [ £||A'£p||^2dr + [ A( —, Qm(t) + || A‘^p||^i)dr 

Jo Jo Jo Co 

< Aoo(t) [ A{Ym{T))dT, for m > 6, (4.6.18) 

Jo 

where (14. 3. 311) has been used in the last inequality. 

Note that 

A*(i) = 4i- A4, A.(”) = Ah _ + ?H|- 4a> 

p p^ pZ P P P^ pA 

Thus, for m > 6, it holds that 

IJ 2 I < TV f [ 7P\2^A‘^p\^dVrdT + Ce^ T || Z“(7p5rpA^(T)) 

16(2/r + A) Jo Js Jo IP 

< pofd 'A [ [ + Aoo(t) [ A( —, Q™(t) + || A>||^i + £|| A>||^2)dr 

16(2/i + A) Jo Js Jo Co 

< x 6(2/!+ A) / _^7Pl-2^“AVpdVrdr + Aoo(t)y A(y^(r))dT, (4.6.19) 
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and 


\J 3 \<e{l\\Z-A‘^pfdry 


|Z“V^p- + ||ZV^pZA‘^(-)f + ||V‘^pZ“A‘^(-)f dr 

P P P 

rt r ^ rt 


}' 


< 


16(2^ 


£7p|Z“A>|2dV.dr + Ce2A(l Q(l) + Q™)^ (\\V^p\\^^. + \\V^pV^vf^. 


+ \\Z<^V^p\\l\\A^v\\l + ||Z“V‘^p||Ll|A>||2 + \\A^{v,p)rn. + \\^^iP,v)\\n^)dr 


< 


16{2p 


< 


/ f jp\Z<^A'^pfdVrdT + CA{-,Q{t) + Qrn) [ e^WvMn-+ HymiT))dT 
^ Jo Js Co Jo 

-^—— [ [ 'yp\Z°‘A'^p\^dVrdr + Aoo{t) [ A{Ym(T))dT. 

P + A) Jo J 5 Jo 


(4.6.20) 


16(2Ai- 

Similarly, for m > 6, it holds that 

\Ja\ < j^^p\Z^A^p\^dVrdT + A^{t)j^ A{Ym{T))dT. (4.6.21) 

Next the expression 

div‘^(p(i9f + V ■ A/'^)v) = gdtdw'^v + QVyA/yddy'^v + -— ^pl4-^3div‘'’t' + V'^gdfv + \7‘^{gvY\7‘^v, 

z 

implies immediately, that for m > 6, 


Thus, 


||pdiv‘^(p(i9f + V ■ V‘^)v)\\'^ 2 dT < Aoo(t) / A{Ym{T))dT. 


ksl < A) Jo y^P\^“^'^P\'^dVrdT + C \\pdiv'^{g{df +v ■V‘^)v)\\l^ 2 dT 


(4.6.22) 


< 


■ f j 7p|Z“A>pdVrdr + Aoo(t) [ A{Ym{T))dT. 

Jo JS Jo 


16(2/r + A) Jo Js 
Substituting (I4.6.17I) - (I4.6.23I1 into (14.6.1611 shows that 

e||AXr)||^2 + / ||A'^p||^2dr < Aoe||AX0)||^2 + Ao / A{Ym{T))dT. 
Jo Jo 

On the other hand, it is easy to obtain that 

||A>||^. < Ao||(A»(0)||?,. +Ao /‘||5t(A»||^.dr 

Jo 

< Ao||(A‘^p)(0)||^i + Ao / ||A‘^p||^ 2 dr + C'5Aoo(t) / A{Ym{T))dT. 

Jo Jo 


(4.6.23) 


(4.6.24) 


(4.6.25) 


Then, combining (14.6.251) and (j4.6.24l) proves (14.6.151) . Therefore, the proof of this lemma is completed.□ 

We now turn to the most difficult part of L°°-estimates: the control of || V'(;(t)||^i,oo and e\\d 2 zv{t)\\'^. 
It follows from ()4.4.1I1 . (I4.4.3|l . ()4.4.6I1 . (I4.6.2|l . ()4.6.4I1 and ()4.6.15ll that, for m > 6 

l|v^(t)||?,.„. < A(l, |M0l«3.=o)(||div^^(t)||^..^ + ||5„(t)||?,.„. + Mt)\\h.^) 

Co 

< WA^pWI. + A(l, Q5(i))(l + WVpfuY + Ao||5„(t)||^c» 

Co 

rt 


< A(y6(0)) + Aot sup A(ym(r))(l+ [ ||Vp(T)||^5dT) + Ao||S'„( 
0<T<t Jo 




(4.6.26) 
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where the following elementary estimate have been used 


Qbit) < 25 ( 0 )+ (7 [ Ye{T)dT, 

Jo 

||Vp(t)||2,4 < ||Vp(0)||?,4 + [ ||Vp(r)||2,.dr. 


(4.6.27) 

(4.6.28) 


Similarly, one can get that 

eWd.Mil < A(-, |Mi)l«3,o.)s(||Vdiv‘^r;(t)||L + ||V^„(t)||L + ||Vz;(t)||^.,.) 

Co 

< Aoll+ A(l, Q,(t) + WVvmDeWVpWl. + Aoe(||5„(i)ll«c=o + IIV5„(t)||L) 

Co 

< A(y6(0))+Aot sup A(F™(r))+Ao£(||5„(t)||?ii,^ + ||V5„(t)||L), (4.6.29) 

0<T<t 


Therefore, it suffices to bound ||S'n(i)||^i,oo and e||VS'„(t)||^. 

However, it seems difficult to bound ||5'n(t)||^i,oo and e||V5'n(t)|l^ by using the approach in [32] 
directly. This is because (Vy, 0)‘(V‘'’div‘^n • N) must appear in (j4.4.12l) for compressible flows, but we 
can not expect to bound uniformly || Vdiv‘^u|||^i,oo and e|| Vdiv‘^ti||^ 2 ,oo■ Thus, we propose to estimate 
Wn, which eliminates V'^div'^n in its equation, instead of S^. Indeed, it follows from (14.5.211 that, for 

TO > 6, 

||^n(i)ll«uo. < Ao|K(t)||?,..„ 

Co 

1 7* 

< A(—Fm(0)) + Ao||a;„(t)||^i,oo + C / A(Ym{T))dT, (4.6.30) 

Co Jo 

and 

£||V5„(t)||^ < Ao£||V a;„||L + A(-, |M0l«3.=.)(e||Vu(t)||?,.,. + s||n(t)||?,.,.) 

Co 

< Aoe|| Vwnll^ + II A'^pIII^i + A(—, Q5(t))(l + e|| Vp||^4) + Ao||a;n||^i.oo 

Co 

1 7* 

< Ao£||Vwn||L + Ao||w„||^i,cx, + A( —,F^(0)) + Aq / A{Y^{T))dT. (4.6.31) 

Co Jo 

Thus, it suffices to bound ||wn|||ii,oo and £|| Va;n(7)|l^. However, the disadvantage of the term Wn is 
that it does not vanish on the boundary. But fortunately, the expression of Wn in (14.5.211 and the dynamic 
boundary condition (11.3111 yield that, on the boundary 

= -2n{(Vy,0)‘at/i- (V‘^N)‘i;^}. 


Thus, we introduce 


c„ := + 2n{(Vy, OYdtV - (V‘^N)‘z;} (4.6.32) 

which vanishes on the boundary. Then, for to > 6, it follows from (j4.6.2p . (14.6.2711 and (j4.6.28ll that 

||^n||?,i,oo < ||Cn||?,C=» + A(l, |Mt)l«3,=o)(l + < IICn||?,C=o + A(l, Q,{t)) 


and 


Co 


Co 


1 7* 

A ||Cn|lw3-“ 4" A( — ,ym(0))+Ao / A{Ym{T))dT, 
Co Jo 


£||Vw„||L < £||VCnIlL + A(-, |/l(t)|^a,=o)£(l + IIv«(i)||^) 

Co 

<e||VCn||^ + A(l,r^(0)) + Ao [ P(YmiT))dT. 
Co Jo 


(4.6.33) 


(4.6.34) 
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As a consequence, it remains to estimate ||Cn(t)||^i,cx, and e||VCn(t)||^- Moreover, by using the 
following lemma, one can see that actually it suffices to derive these i°“-estimates near the boundary, 
that is ||xCn(t)||^i,oo and e||xVCn(t)||TC with x compact supported and equal to 1 in the vicinity of z = 0. 


Lemma 4.25 For any smooth cut-ojf function x such that x = 0 in a vicinity of z = 0, it holds that for 
m > k 2, 

WxfWn^.^ < (4.6.35) 

This Lemma follows from the Sobolev embedding and the fact that the standard Sobolev norms are 
equivalent to the conormal ones away from the boundary. In fact, Lemma 14.251 and (14.6.21) imply, for 
m > 6, that 

IICn||?ic=. < ||xCn||?ii,=. + 1^11?,.,=. < ||xCn||?,c=o + A(l, ^^(O)) + f A{-,Y^{r))dT. (4.6.36) 

Co Jo Co 


One can obtains the following lemma by the same argument in [32] but using ||xCn||^i oo instead of 

WxSAl^,^-- 

Lemma 4.26 For m > 6, it holds that 

||Vu(t)||^i,oo < A(—,ym(0))+t sup A( —,ym(T)) . (l+£^|Z'"h(t)|i + f ||Vp(T)||^ 5 dry (4.6.37) 

Proof. As indicated by (I4.6.26|) . (14.6.301) . (14.6.331) and (14.6.361) . we need only to estimate ||xCn||^i,oo. 
The main difficulty of this estimate is to handle the commutator between Zi and We shall use a 
normal geodesic coordinate system in the vicinity of the boundary which gives a simpler expression of 
Laplacian|32]. Here, we can use this coordinate system because L°°-estimate does not require the highest 
regularity of the boundary. Let us define a new parametrization at the vicinity of the boundary 

4'(t, •) : S X (—00,0) —)► D{t), 

(y, z) ^ (y, h{t, y)) + zn^{t, y), (4.6.38) 


where D{t) is defined in (11.251) and n*’ is the unit exterior normal rf’{t,y) = (—9i/i, —92/^, 1)/|N|. Note 


that 


/ 1-Z9i(^) -Zd2i^\) 


D^{h-) = 


-Zd^{^\) 1-Za2(^|) 


IN 

d^.h I 


|N| 

a^.h 


(4.6.39) 


yai/l + 2;ai(i^) 92/1 + 2^2 (j^l) 1^) 

which is of the form Mq + R with |i?|oo z\h\ 2 ,oo and 


Mo = 


/ 1 
0 


n dih 

6 

1 —§2h. 

^ INI 


y dih d2h 


is invertible. This yields that 4'(t,.) is a diffeomorphism from x (—(5,0) to a vicinity of d^t for some 
(5 > 0 which depends only on cq > 0. By this parametrization, the scalar product in induces a 
Riemman metric which is given by 


9{y.z) 


^y(y,2) 0^ 


and the Laplacian in this coordinate system is of the form: 


Ag/ = d^^f + ia^(ln \g\)dj + Ag/, 


where |y| denotes the determinant of the matrix g and Ag/ is given by 

l<^.g<2 
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which involves only the tangential derivatives. 

In order to use this normal geodesic coordinate system, one first localizes the equation for (n- Set 

= x(z)a; (4.6.40) 

where x(.^) is given as xi^) = k( ) with k{z) be a smooth function of compact support such that 
k{z) G [0,1] taking the value 1 in the vicinity of z = 0, and (5(co) to be determined later. Note that this 
choice implies that \x’^{z)\ < Aq. 

Applying V^x to yields that in S, 

gOfuj + gv ■ — g.eA'^u} = (4.6.41) 

where 

X dfv — V^g X {{v ■ V‘^)v) + guj ■ — gdiv'^vuj. (4.6.42) 

Furthermore, it follows from (14.6.401) and (14.6.411) that 

gdfuj^ + gv ■ - neA^uj^ = F^x (4.6.43) 

where := F^ + xFu) with 


Ff = g{Vzdzx)^ — — /xeA'^xw. 


(4.6.44) 


Since F^ is supported away from the boundary, it follows from (13.221) , p.23l) and (14.6.351) that 

\\Fnl^.^ < A(-, + \\{v,p)ml.^,^)\Ht)\\l.. < A(-, Q,{t)). (4.6.45) 

Co Co 

Next, we define uS^ := w^(t, •) o vk). This change of variable is well-defined if we choose (5(co) small 

enough such that Dd! is invertible. Then, solves in S the convection diffusion equations 

gdtuj'^ + £ib • Vw'*' - ne{dzz + ^9^(ln \g\)dz)uj^ = F^x{t, o vk) -|- g^eAgU;'^, (4.6.46) 

where 

b = {Dd))-\v{t, o 4-) - (4.6.47) 

Set 

2/: z) = n^(t, y)iuj'^(t, y, z) X N^(t, y)), (4.6.48) 

with = {—dih,—d 2 h,l) and 11^ = Id — n*” (g) n^. Note that 11*’ and n*’ are independent of 2 . This 
yields that solves 


where F^ is given by 
with 


gdtcjn + 6’b • Wujn - pe{dzz + ^^^(In \g\)dz)ujn = 


^ _ 171 ^ 

(jJ •> 


FJ = n^’F^.x X N" + FJ’^ + FJ’^ 


Fj’i = g{dtU^ + hy ■ Vn*')^'’' X N*’ + yn{w'^ • {dt + b,, • Vy)N*’}, 
F^'^ = -iieIl\A~giv'^ X N*’). 

Define 

= ixv)(t,<^~^ o 4'). 

Note that V'^N*' = VN** is independent of 2 . Thus n*'{(V‘^N*’)*u'^} solves 


(4.6.49) 

(4.6.50) 

(4.6.51) 

(4.6.52) 

(4.6.53) 


yatn*’{(V‘"N'’)‘r;'*'} + yb • Vtf{(V‘^N^)‘z;'^} - pe{dzz + \dz{ln |y|)a,)n*’{(V‘"N*’)‘i;'''} = F„'^, (4.6.54) 


n*'{(V‘"N*')‘(xF„ + F})} + + F^’^ 


where 


(4.6.55) 
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with 


= -V> + (m + A)eV'^div‘^?;, 

F} = QiyzdzX)v — ~ 

F^’^ = gidtU’^ + hy ■ Vtf)(V‘^N'')^'^ + 0n{(9t + hy ■ Vy)(V‘^N)‘z;'^} 
F^’^ = -nsIl\iV^N^YA~gv'*). 


(4.6.56) 

(4.6.57) 

(4.6.58) 

(4.6.59) 


Set 


C = - 2tf {(V'^N'’)*?;'^} + 2n'’((V,„ OYdtv). 

It follows from (I4.6.49|) and (14.6.541) that Cn solves in S 


(4.6.60) 


QdtCn + Pb • VCn - tJ.e{dzz + l5l)52)Cn = -Pj + 


(4.6.61) 


where 


= gdt{Il\Vy, oydtv'^) + gh ■ V(n''(Vj„ oyotv'^) - gieidzz + -^.(In |g|)a,)(n'>(Vy, Oydtv'^), (4.6.62) 
with defined by x^(i; o '!')■ In order to eliminate the term ii9z(ln , we define 


r'^ — _1_A^ 

I I 1 • 

Ish 


Plugging into (I4.6.61L one obtains that 


where F^' 


gdtC + • VC - g.edzzC = + F^) =: F, 


Py =-{gdt{\g\ «) + (?-bV(| 5 | ^ - nedzzi\g\ ") - ^ )Cn 


\9\- 


Moreover, the boundary condition reads as 


C = C = Cn = 0 . 


(4.6.63) 


(4.6.64) 


(4.6.65) 


Since Sobolev conormal spaces are invariant by diffeomorphisms which preserve the boundary (see Lemma 
9.5 in [35]) and |n — n^| + |N — N^| = 0{z) in the vicinity of the boundary, one obtains that 

||Cn||?,i,=o < Ao(||Cll«U<» + + llulll,.,.) < Ao||Cll«c=o + A( 1 ,Q 5 ( 0 ) 

Co 

< Ao||Cll«i,-+A(P™(0))+ / A{Ym{T))dT (4.6.66) 

Jo 

which implies that it suffices to estimate ||Cn IIhi.=o- Now, we will estimate ||Cnll|ii,oo by taking advantage 
of (14.6.631) and (I4.6.65|) . First, since 


63(t,y,0) ^^{D'^)3yivj -9t4'j)(t,2/,0) = 0, 
i=i 


then (14.6.6311 can be rewritten as 

0{t, y, 0)(9t + bi{t, y, 0)5i + b 2 it, y, 0)^2 + zdzb^it, y, 0))C - y^dzzC^ 
= S - [g{t, y, z) - g(t, y, 0)]9tCn “ 2^’ 0)]^*Cn 


(4.6.67) 


- Q{t,y,0)[b{t,y,z) - zdzb3{t,y,0)]dzCn “ [Q{t,y,z) - Q{t,y,0)]zdzb3{t,y,0)dzCn =■ G. 


(4.6.68) 
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It follows from Lemma 16.21 in the Appendix that, for m > 6, 


IlCn < IIC(0)||«n=. + ||p-ioo||G||„.,=.dr 


Al + ||p-i|U)(l + ^||Z,(e,6,a,63)||L)|ieil«^-rfr (4.6.69) 

•^0 i=o 


< 


11^(0)11 


HI,. 


■A, 


: TllGIlHno 

Jo 


.dr + / A{ — ,Ym{T))dT. 

Jo Co 


It remains to estimate ||G||Hi.oo(iT. First, it follows from (j4.6.42|) . (14.6.451) and (I4.6.50I) - (I4.6.52I) that 

< A(l, Q 5 (i))((l + ||Vp||^,,^)||Vr;||?in=. + ||+ e^HVz;||^ 3 ,=.) 

Co 

< + e^\\V^v\\n4^v\\n^), m > 6, (4.6.70) 

Co 


where Lemma [4. 2 31 has been used several times. Next, in a similar way, one can get from (14. 6. 55^ - 114. 6. 59|) 
and Lemma 14.231 that 

< A( 1 ,Q 5 ( 0)[1 + \\yp\\H^,^ + ||Vr;||^,,^ +£ 2 ||Vdiv%||^,.^ +e 2 ||v«||^ 3 ,=o] 

Co 

< e^llA‘^ p||^2 + A( —, QG{t) + ||A‘^p||^i + ||Vti||^i,oo)(£:^||Vp||^5 + e^||Vu||^3,oo) 

Co 

< Ai-,Y^{t))il + e^\\v4\\n4^v\\n4, rn>6. (4.6.71) 

Co 

It follows from (14. 6. 621) . (14. 6. 641) and Lemma [4.231 that 


+ \\\g\iF^rn.,^ < A(f ,r™(t))(l + e^\Z^h{t)\l), m > 6. 


Co 


Therefore, combining (j4.6.7ni) - ()4.6.72l) leads to that 

< A(l,y™(t))(I + e^\Z"^h{t)\l + e4vMn4^v\\n4, rn>6. 
Co 2 


(4.6.72) 


(4.6.73) 


Next, using the Taylor’s formula and the fact that () is compactly supported in z, one obtains that 
II {Q{t, y, z) - g{t, y, 0))5tCn IIhu- 


< AollClI^uoo + \\Z{g4y,z) - p(t, y, 0))||^||a*C||^ + \Mt,y,z) - p(f, y, 0))Za*C|| 


a^|[ 2 
I oo' 


< AollClI^uo. + ||2y|li,|lCll^ + \\d.g\UY^^2dtC\\ 

Due to (1^ . the following inequality holds for |a| > 0 

llT^^^eilL < Iiv(^z“e)||i||^z“eil2 < l|2“eil^ 

I — z 1 — z 1 — z 

Therefore, plugging (14.6.751) with |q;| = 2 into ()4.6.74|) . one gets from (14.6.81) that 


I OO 


(4.6.74) 


(4.6.75) 


|l(y(t,y,z) - git,y,0))dtCn\\W’'=° < H — ,Ymit)), m > 6. 

Co 


(4.6.76) 


E 

2 = 1,2 


||y(&,(t,y,z) - 6*(t,y,0))a*Cn IIhi.oo < A( —,y,„(t)), m > 6. 

Co 


Similarly, 


(4.6.77) 
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Since b^{t,y,Q) = 0, it holds that 

\\Q{t, y, 0)(b3(t, y, z) - zd^b^^t, y, 0))d^Cn llwi.- 

<\\9\\w^M\\^b3\\lKni^.^ + mh^ 

+ \\{b3{t,y,z) - zd^b3{t,y,0))Zd^C^\\l^) (4.6.78) 

<A(l,Q5(0)(l + l|Vn||?,.,^)||Cll|?,c=.+A(l,Q5(i))l|az.&3|lLll7T^^^^C^IlL 

Cq Cq yi z ) 

<A(-,Q5(0)(l + l|Vn||?,.,<.)||e|||,.,<. 

Co 

< A( —,47„(t)) + A(—,47„(t))|19^;,53||^||Cn 11^4, for m > 6. 

Co Co 

It follows from 1)4. 6. 391) and (14.6.471) that 

||^^..63||L<l|5..^-N|lL + A(l,|h|?,4,o.) 

Co 

< AoP^div'^'ull^ + Ao||V?;(t)||^ + A( —, |h|^4,oo) < A( —,F^(t)), m > 6, (4.6.79) 

Co Co 

where one has used that 


dzz{v ■ N) = dz{dzV ■ N) + dz{v ■ 5^N) 

= dzidzifidiv'^v - divi - 52 ^ 2 )) + dz{v ■ dzN) 


Thus, 


\\g{t,y,0){b3{t,y,z) - zdzb3{t,y,0))dzCn\\n^.^ 


< A{—,Yjn{t)), m>6. 
Co 


(4.6.80) 


Similarly, note that {g{t, y, z) — g{t, y, 0))2; = 0{z‘^) in a vicinity of the boundary, one gets that 


\\iQit,y,z) - g{t,y,0))zdzb3{t,y,0)dz:C^\\l^i,^ < A( —,47„(t)), m > 6. (4.6.81) 

Co 

Thanks to (|4.6.73|) . (I4.6.76L (I4.6.77L (14.6.801) and (14.6.811) . one has that 

||G||?,4.„ < A(l,y™(t))(l+£2|z-h(t)|2 +£2 ||v\||^,), m>6. (4.6.82) 

Co 2 

Substituting (14.6.821) into (14.6.691) and using (14.3.311) . one gets that 

llC^WII«u=»<|ie(0)||?,uoo+ f A{-,Y^m^+e^\Z^h{t)\l)dr, m > 6. (4.6.83) 

Jo Co 2 

Then, (14.6.371) follows from (14.6.261) . (I4.6.30p . (14.6.331) . (I4.6.36L (14.6.661) and (14. 6. 831) . Thus, the proof of 
this lemma is completed. □ 

Lemma 4.27 For m > 6, it holds that 


£\\9zzv{t)\\l^ < t sup A —,r^(r)+£2|Z™/i(T)|iV (1 + / ||Vp(T)||^5dT) 

0 <T<t \Co V Jo 

+ A(1,F™(0)). (4.6.84) 

Co 

Proof. We shall reduce the problem to the estimate of ejl^jCn Iloo- First, it follows from (14.6.291) . (j4.6.3ip . 
(I4.6.34P and (14.6.371) that 

< t^sup^A(^,F^(T))(^l+ £2|Z'"/i(t)|| + J ||Vp(T)||^5dr) 

+ A(l,y„(0))+Ao£||VC„||L,"i>6. 

Co 


(4.6.85) 
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Note that |n — n^| + |N — N*”! = 0(z) in the vicinity of the boundary, one gets from the definitions of 
Cn and Cn that 

eWd^Uml < Ao(£||9.Cn IlL + IlK 

< A( —,rm(0)) + Aot sup A{—,Ym{T))(l + e'^\Z"^h{t)\\ + [ ||Vp(T)||^ 5 dT') 

Co 0 <T<t Co V ^ Jo ' 

+ AoellS^Cn ||oo,m-> 6. (4.6.86) 


It remains to estimate e||9^Cn IlL- Since Q solves (14.6.631) in S with the homogeneous Dirichlet boundary 
condition, one can use the one-dimensional heat kernel of z > 0 


G(t,z,z') 




(z - z’f 
4^et 


) - exp( 


(z+^ y 

4/xet / 


to obtain 


+ 00 


/ d^G{t,z,z)Q({s)dz 

(•t /* + oo 


^/ed^G{t - T, z, z')(5' - pb • VC^ - (p - l)dtC,Z)dz'dr. 


0 ^0 


It follows from (0) = 0 and integration by parts that 


\/e||5^C^lloo < \/e||9^Cn (0)llc 

It then follows from (14.6.731) that 


1 


y/t^ 


'o y/t-T 

r.t 


ll-Sllwi.oodr 


< 


^ A(l,y„(t)) fl + e\Z”^h{t)W + ^e\\V\d,V~e\\yv\\ 


y/r^ 


T Co 




dr 


<tA (-^Y^{t) + £^\Z’^h{t)\\ 

\co = 

where one has used (14.3.311) in the last inequality. Since 63 (t, y, 0) = 0, so 

1 

Co 

Finally, it is easy to obtain that 


pb. vciiy < Ao(ii6|iy + iiv6|iy)iiciiy4,. < A(-,F^(t)). 


Co 


Substituting (I4.6.89l) - (|4.6.91l) into (14.6.881) . one gets that 


(4.6.87) 


+ lipb. VC^IU + ||(^?- l)5tC^||oo). (4.6.88) 


(4.6.89) 


(4.6.90) 


(4.6.91) 


e\\U^\\l<V^m^m\lo+tAi-,Y^it)+e^\Z^Ht)\l). (4.6.92) 

Co 2 

Now (|4.6.84l) follows from (I4.6.92L (14.6.851) and (14.6.861) . Therefore, the proof of this lemma is completed.□ 


4.7 Proof of Theorem 14.11 

Based on the estimates obtained so far, we can complete the proof of Theorem l4.1l in this subsection.First, 
(14.0.11) follows easily from conservation of mass. 
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In order to prove (I4.0.2I1 . one chooses ^ > 0 large so that 

dz^{0,y,z) > 1, 


(4.7.1) 


where A depends only on the initial data |h(0)|^f3. For two parameters i? > 0 and cq to be chosen 
0 < — i?, define 

Co ’ 


]; Om{T) < R, d^ip{t) > Co, |h|«3,=o + |Vyh| [^]+i < —, 

Co 


- arp|..o i f. 0 < Jjf S e(‘) S ■>Co,w e [o,r|}. 


(4.7.2) 


The restriction T < guarantees that the fluid moves a distance no more than idg. Then flo.i, 

fifc i, fc = 1, • • • , n is still a valid covering of Q,t with t <T. 

?^ow define 

=Afm{p,v,h){t) 

sup |l + ||(F'",Q'")(r)f + |(h,v^V,h)(T)|^„ + ||(Vp,V^;)(T 

0<T<t '■ 


l^m-2 


+ ||A‘^p(r)||^. + ||V(p,^)(r)||^.„. +s||(divVV>5„)(r)||^,„_. +e|Z-h|i +e||V^^(r)||^., 


+ £||V(p,z;)(r)||^„._i +e||AXT)||^2| ||Vp(T)||^„_idT + ||V?;(T)||^™_idr) 

+ e f ||Vn(T)||^™dr + e / \\V^v{T)\\l^m-. 2 dT + f || V\(t)|1^,„_i dr, 

^0 -'0 Jo 


(4.7.3) 


/o -/o 

which implies immediately that 

Ym{t) + M{t) < AoA(A/m(t)), A/)„(t) < AoA{Om{t)) and Om{t) < AoA{Afm(t)). (4.7.4) 

It follows from Propositions 14.17114.161 that 

IKF'", Q'")(t)f + \{h, V^Vyh){t)\l^r. + IKdiv^n, V>, SnmWn^-^ + £|l(div‘"^^, V>, 5„)(t)||?,„_, 

+ £|Z™h||+^ ||Vp||^,„_idr + \\Vv\\^m-idTy +e \\Vv{t)\\ 1^„ + \\V^v{T)\\'^m- 2 dT 


+ s^ 


■ [ \\VMr)\\n--idT < AoA(Ar„(0) + ||Vn(0)||^,„_3) + tiAoA{Afm{t)). 
Jo 

Then, as a consequence of (j4.7.5p . Lemmas l4.24[ 14.251 and [T??fl one obtains that 


Therefore, 

Moreover, 


em(0< A(Co,0„(O) + ||V^;(O)||^,„_3)+t5A(Co,i?), Vte [0,T^]. 


Cc 


exp(rA(i7)) < g{t) < Co exp{TA{R))yt € [0,r], 


< |h(0)||,3,o. + |v,h(0)|^[^,+, +TA(i?), vt e [o,r], 




dMt) > 1 - / \\dt^vit)\\L^dT > 1 - TA(i7), Vt e [0,T], 

Jo 

once the choice of A > 0 so that ((4.7.111 is satisfied. Finally, note that 

-dtp{t) > -d'^po - [ ||5t(9|’p(T)||dr > -d'^po - TA{R), Vt e [0,r]. 
Jo 


(4.7.5) 

(4.7.6) 

(4.7.7) 

(4.7.8) 

(4.7.9) 

(4.7.10) 


(4.7.11) 
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Taking R = 2A(C'o, (9^(0) + ||V?;(0)||^m_i), in view of (I4.7.7|) - (I4.7.11I1 . one gets that there exists T» > 0 
depending only on Cq, do and the initial data Omifi) + ||Vz;(0)||^m-i (hence independent of e and a) 
such that for T < niin(r*,T'^) 

'e^{t) < 2A(Co,0™(O) + ||Vu(0)||^_0, 

< d,ip{t) > 2co, |/i(t)|«3,=o + |Vy/i(t)|^[m]+i < ^, vt e [0,r]. (4.7.12) 

> ^, 0 < 5 ^ < g{t) < 2Co, 

Of course, it holds that T, < Tf. Then taking Ta = T,, one obtains (14.0.211 and closes the a priori 
assumption 114.0.41) . Therefore, the proof Theorem HTI is completed. □ 


5 Proof of the Main Results 


Proof of Theorem 11.51 We shall combine the a priori estimates obtained to prove the uniform existence 
result. Fix TO > 6, consider initial data {pq,Vq, h^) € such that 

Zm(0):= sup + ||Vug||^„-i) < Co, 

eG(0,l],crG[0,l] ^ ^ 


and 

0 < ^ < pg < Co, (5.1) 

For such initial data, we are not aware of a local well-posedness result for CH) and dEll) and (HJ). 
Thus, we shall prove the local existence by using energy estimates and a classical iteration scheme. 
Since (pq, Uq, Fq) € there exists a sequence of smooth approximate initial data (pg’^, Ug’^, Fg’^) G 

^N^ap(^ being a regularization parameter), which have enough spatial regularity so that the time deriva¬ 
tives at the initial time can be defined by Navier-Stokes equations and the boundary compatibility con¬ 
ditions are satisfied. We construct approximate solutions, inductively, as follows 

(1) Define = Ug^, and 

(2) Assuming that was defined for k > 1, let (p^,u^) be the unique solution to the following 

linearized initial boundary value problem: 


Pf + div(p'^u^ ^) = 0 in (0,T) x fit, 

p^Ut + p^u^~^ ■ Vu^ + Vp^ = eAu^ + eVdivu^ in (0, T) x fit, 


dtF^ + ■ VF^ 

{p^,u^ icfc 

with boundary conditions 


= 0, on dfit, 
0 


F>^)\t^o = {pr,<’'',Fon, 


(HH). 


with 3^ < pI’^ < |Co, 


(5.2) 


Since p^, and F^ are decoupled, the existence of the global unique smooth solution (p^,it^) of 
(lOll with 0 < p^ < oo can be obtained by using known results. For example, one can first obtain a 
local existence result for (15.21) by using local existence result in [JD] without surface tension and [531IM] 
with surface tension. Since dSJ) is a linear problem, one can extend the local solution to a global one. 
Although the local existence result of [40] without surface tension and |63l |64] with surface tension is 
proved in Sobolev-Slobodetskii spaces W*’ 2 (see m) with I > 0, but it is enough for us to proceed our a 
priori estimates in Section 4 by taking I = 3to. 

Similar to Section 1, let {p^, v^, h^) be the localized version {p^, , F^). Applying the a priori estimates 

given in Theorem 14.11 and by an induction argument, we obtain a uniform time Ti > 0 and constant 
Ca > 0(independent of e,cr and 5), such that it holds for {p^,v^,h^), k>l that 

0™(p^u^h'=)(^) < 03 , Vf e [0,ri], (5.3) 


and 

< p"(i) < 200, vte[ 0 ,ri], (5.4) 

ZOo 

where Ti and 03 depend only on 0o, 0o and do. Based on the above uniform estimates for {p^,u^,F^), 
using Lagrangian coordinates and considering the difference of — p^, , 0^+^ — F^), one can 
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prove, by using the energy methods, that there exists an uniform time T 2 (< Ti)(independent of e,a and 
S) such that converges to a limit as k ^ +oo in the following strong sense: 

( (p^u'=) ^ in L°°{ 0 ,T 2 ;L^), 

< ^ in L 2 (o,T 2 ,L 2 )^ 

[ F’^(t, xi; ■)) ^ in L^iU). 

Then, it is easy to check , F’^’^) is a weak solution to the problem (ll.lll . (11.511 and (HSI with 

initial data {p^’^ ,Uq’^, F^’^). Then, by virtue of the lower semi-continuity of norms, one can deduce from 
the uniform bounds (15.31) and (15.41) that ,u^’^ , F'^’^) satisfies the following regularity estimates 

< Cs, vt e [0,T2], (5.5) 

and 

^ < 2Co, Vte[0,r2], (5.6) 

zOo 

Based on the uniform estimates (ESI and (ESI for (p'^>'^, F®’^), we pass the limit (5 —>■ 0 to get a 

strong solution (p'^,u®,F'^) to (11.11) . (I1.5|l and (11.611 with initial data {pq,Uq,Fq) satisfying (|5.1I) by using 
a strong compactness arguments. Since the compactness arguments are almost the same as the ones need 
for the proof of Theorem II. 101 below, we shall not give more details here. Moreover, the uniform bounds 
(11.471) and (|1.48|) are immediately results of the lower semi-continuity of norms. Changing the variable 
into Lagrangian coordinate and using energy method, it is easy to prove the uniqueness of {p^,u^ , F^) 
since we work on functions with enough regularity. Taking Tq = T 2 and Ci = C 3 , we complete the proof 
of Theorem [TE □ 


Proof of Theorem 11.101 For any fixed cr > 0, it follows from Theorem 11.51 for any e £ (0,1], that 


< Cl, Vf G [0 ,Fq], 


(5.7) 


and 

i^<g%x,t)<2Co: VtG[0,To]. (5.8) 

zOo 

Thus, {p^,v^) is uniform bounded in F°°([0, T]; F™), V(p'^,u^) is uniform bounded in F°“([0, F]; F™“^), 
h® is uniform bounded in F°°( 0 ,T; F™), dt{p‘^,v‘^) is uniform bounded in F°°([ 0 , T]; F^) and dth’^ is 
uniform bounded in F°°(0,T; L^). Then, it follows from the strong compactness arguments(see [12]) that 
{p^,v'^) is compact in C{[0,T]-, and h® is compact in (^([O, F]; F™”^). In particular, there exist a 

sequence — >■ 0-|- and (p, u) G (^([O, F]; F™“^) and h G C{[0,T]; F[^~^) such that 


{p^-,v^-)^{p,v) in C'([0,F];F™-2), 
h’^^h in C'([0,F];F™-i), 


cLs Sfi —y 0“h, 


(5.9) 


Denote (p, u, F) the global version of (p, v, h). Applying the lower semi-continuity of norms to the bounds 
Eli) and EH), we obtain the uniform estimates (11.501) and (11.511) . Then, based on (15.91) and the uniform 
estimates (I1.50L (11.511) . it is easy to check that (p, u, F) is a solution of the Euler equations with the free 
surface 


dtp + div(pu) = 0 , 


X G fit, t G (0, Fq] 


]pdtu + p(u ■ V)it -t- Vp = 0, 
with the boundary conditions 

dtF + u ■ VF = 0, and p = Pe — crF, x G dflt, 
where fit is the domain occupied by the fluid on time t > 0 and the boundary of fit is given by 

dflt = {x G F{x,t) = 0}. 


(5.10) 


(5.11) 
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By using (I1.50L (I5.9|) and the anisotropic Sobolev inequality (14.1.311 . we have that 

sup - ■!;)|||oo < sup (||V(/" - 0. 

te[o,T] te[o,T] ^ ' 

and 

\h^ — h\y/l,ao < \h^ — ^ 0, 

as En —t 0+. Hence, (ll.52|) is proved. 

Consider {p^,u^,F^), i = 1,2 two solutions of (I5.10|) . (15.111) defined on [0, T] with the same initial data 
and the regularity stated in (ll.50|) and (11.5111 . We shall prove that (p^,u^,F^) = (p^,u^,F^). Applying 
the standard energy method to the difference (p^,u^,F^) — (p^, , F^) in the Lagrangian coordinates and 

using the uniform bounds in (ll.SOp and (11.511) . Therefore, {p,u,F) is a unique solution to Euler system 
with free surface (I5.1()l) - (|5.11ll . And the uniqueness yields immediately that the whole family (p®, u®, F’^) 
converges to {p,u,F). Thus, the proof Theorem (TTU] is completed. □ 

Proof of Theorem 11.111 Based on the uniform in both e and a estimates in Theorem 11.51 one can 
prove Theorem 1 1.11 1 bv using the similar strong compactness argument in Theorem ll.lOl The details are 
omitted here for brevity. □ 


6 Appendix 

Let S{t,T) be the (7° evolution operator generated by the following equation 

[dth + hi{t,y)dyih + h 2 {t,y)dy 2 h +zh 3 {t,y)dzh]-ed{t,y)dzzh = 0, z > 0, t > t, ( 6 . 1 ) 

with the boundary condition h{t,y,0) = 0 and with the initial condition h{T,y,z) = hQ{y,z). The 
coefficients are smooth and d{t, y) satisfies 

Cl < d{t,y) < — (6.2) 

Cl 

for some positive constant ci > 0 . 

Then we have the following Lemnas whose proof can be found in |53) . the details are omitted here. 
Lemma 6.1 It holds that, for t > t > 0 

l|5'(t,r)/io||L“ < ll^olU”, (6.3) 

\\zd^S{t,T)ho\\Lo^ < C{\\ho\\L'^ + \\zdz:ho\\L’^), (6.4) 

where C > D is a uniform constant independent of the bound ofd(t,y) and bj{t,y), j = 1,2,3. 

Lemma 6.2 Consider h a smooth solution of 

a{t,y)[dth + bi{t,y)dyih + b 2 it,y)dy 2 h +zb 3 {t,y)dz:h]-edz:zh = G, z > 0, h{t,y,0) = 0, (6.5) 

for some smooth function a(t,y) satisfies (16.211 and vector fields b = ( 6 i, 62 , 63 )*(t, p). Assume that h and 
G are compactly supported in z. Then, one has the estimate: 

\\h\\n^.-<\\ho\\n^.^+ T |li|U»||G||„c=odr 

Jo a 

+ Al + I|-|1l^)(1 + E b)\\l^)\\h\\nr^dr. (6.6) 

Jo a ^ 
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